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I. Collective Oscillations and Characteristic Electron-Energy Losses 


By D. Gapor 


Imperial College, London* 
[Received September 1, 1955] 


§ 1. THE NaTURE OF THE COLLISION PRocESS 


SHARPLY defined ‘characteristic’ energy losses of medium-energy 
electrons in thin films, sometimes approaching the appearance of 
line-spectra, have been observed by a great number of authors.t No 
acceptable interpretation was forthcoming until Bohm and Pines (1952) 
and Pines (1954) suggested that they might represent the transfer of one 
or several vibrational quanta from the fast electron to the collective 
oscillations of the electron gas in metals.f 

Since that time two pieces of experimental evidence have been found 
which appear to give strong support to the ‘ collective’ interpretation. 
One is the recent observation by Watanabe? of a strong correlation between 
energy loss and direction of the electrons, which, as will be discussed below, 
appears to give a direct verification of the dispersion law w(k) of the plasma 
waves. The other is the observation, by Gabor and Jull (1955) of the 
strongly non-linear increase of the collision cross section with film thickness. 
Neither of these is compatible with an interpretation by atomic collisions, 
but further evidence is required before rejecting other alternatives. 

It will be attempted to give a simple but adequate theory of the 
interaction of fast electrons with thin films of solids, sufficient for 
covering past experiments and for suggesting new, crucial tests. 


* Communicated by the Author. 
+ There is no need for exhaustive references as two very comprehensive 


reviews will appear shortly, one on the experimental side, the other on the 
theory. The tirst is by L. Marton, L. B. Leder and H. Mendlowitz, 1955, 
Advances in Electronics (Academic Press), the second by D. Pines, 1955, 
Advances in Solid State Physics (Academic Press). I thank these authors for 
letting me read their manuscripts. 

+ To be published: Quoted by Pines /oc. cit. 
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Interaction of electrons with collective oscillations is a rather anomalous 
phenomenon, because it is a coherent but inelastic collision process. It is 
analogous to the reflection of light by moving mirrors, or by sound waves, 
or to Bragg-diffraction on a moving lattice. The analogue of the 
Cerenkov-radiation, used by Bohm and Pines, though illuminating, 
misses the point of coherence. In fact the wave-character of the electron 
is of no importance in an extended medium, as treated by these authors, 
but becomes so in the thin films as used in all experiments. Moreover, as 
will be shown below, it is not sufficient to replace the electron by a plane, 
monochromatic wave, but a more detailed representation of the electron 
beam is required, which takes account of its chromatic and geometrical 
coherence properties. 

As regards the solid film, we will use the simplest model which appears 
likely to represent its relevant features. If a is the thickness of the film, 
we imagine it built up of rectangular ‘ grains * with transverse dimensions 
b,c. We postulate that each grain is capable of vibrational modes corre- 
sponding to a dispersion law w(/), and that all grains oscillate independently 
of one another. Most of the theory below will be applicable to any dispersion 
law w(k) of longitudinal electric waves, and it is irrelevant whether this is 
derived from the side of the free electron gas, as in the theory of Pines 
and Bohm, or by an approach from the optical dispersion theory of solids, 
as in the recent papers by Hubbard (1955) and by Frohlich and Pelzer 
(1955). It is essential only that there must be a strong dispersion, that is 
to say the wave number k must vary widely for a small change of the 
frequency w. In some places, though, we shall make use of the dispersion 
of the free electron gas, as an illustration. 

The main result of the theory is that it gives two new crucial tests of the 
hypothesis of collective oscillations. One is, that the characteristic losses 
must increase very strongly if the energy-spread of the bombarding 
electron beam is reduced. The other is that, other things being equal, 
the effect must become less if the geometrical coherence of the beam is 
increased. These are striking tests, because they mean that chromatic 
and geometrical coherence of the electron beam have opposite effects. 
Neither would have any influence if the phenomenon were an atomic 
collision process. 


§ 2. REPRESENTING THE ELECTRON BEAM 


While in atomic collision processes it is sufficient to represent the 
incident electron as a plane wave, this is not so in the present case, as the 
bombarded object 6, ¢ need not be small compared with the itera width 
over which the electron beam can be considered as coherent. In order to 
avoid complicated calculations, involving partial coherence, we represent 
the electron beam at the object as the juxtaposition of beamlets. each 
coherent in itself, incoherent with the others. 

The lateral coherence width is, nearly enough, one-half of the de Brogli 
wavelength 4,, divided by the angle which the image of the source (che 
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hot cathode), subtends at the foil. For 10 key electrons .. is O-11A and 
the collimation used by most experimenters is in the limits 10-3-10-4 radian 
hence the coherence widths are of the order 50-5004. This happens to be 
a very important range, and more attention will have to be paid to the 
coherence width in the future. 

Second, we must pay attention also to the ‘ chromatic coherence ’ of the 
beam, that is to say its energy spread. We will use below perturbation 
theory for calculating the interaction with the collective oscillators in the 
solid, and for this we require an interaction time. It would be of course 
quite wrong to identify this with the transit time of the electron through 
the foil; the particle has no place in a theory which deals with wave- 
features of the electron. It will be shown later that the transit time of 
a 10 kev electron through a 604 foil of Al corresponds only to a half-cycle 
of the collective oscillation. One could never account for the observation 
that the energy spread in the line corresponding to a loss of 14:7 ev is 
not appreciably more than the primary spread, of a fraction of a volt only, 
if the electron had ‘ seen’ only half a cycle of the oscillation. It is not 
the thickness of the foil, but the length of the wavetrain that counts. 

In order to simplify the calculations, we will replace the inhomogeneous 
beam by pulses of monochromatic waves, of duration 7. It is evident 
from the uncertainty principle that for a hot cathode of temperature 7’ 
this must be of the order #/k7, but the numerical factor is of great 
importance because the collision probability will be found directly 
proportional to 7. In the Appendix a careful estimation of the best fit 
has been carried out, and the result is 


t=ThlkT. (tT ees aa € 9: 


It is useful to express this in terms of the wavetrain length, tv, where v is 
the velocity of the primary electrons. This is 
ru=Th/kTv=(nh/{ke)B/T =7-28/T cm. 

For a tungsten cathode, as used by all experimenters, 7'=2500°c, and for 
10 kev electrons, B=0-197, this is 5700 angstroms ; many times more 
than the foil thicknesses which are at most a few hundred angstroms. 
In aluminium for instance this wave ‘ sees ’ 34 full cycles, irrespective of 
the foil thickness. 


§ 3. THE MopEL oF THE FoIL 


The Bohm—Pines theory of collective electron interaction has not yet 
been extended to non-uniform structures, and it is difficult to say what 
gaps or structural differences must exist between two , grains > to make 
them independent. The experiments of Jull and the author suggest that 
it is sufficient to introduce a carbonaceous * contamination ’ layer of the 
order 10-204 between two gold layers, or a gold film of similar thickness 
between two aluminium layers for effective separation. It is not our 
purpose at present to go into this question ; it will be assumed that the 
grain sizes are given, and we will discuss the consequences. 


B2 
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A rectangular grain with the boundaries z=0, a, y=0,b, z=0,¢ is 
completely decoupled from its neighbours if the electric field is permanently 
zero at the boundaries. Thus the potential field ¢ can be represented by 
a sum of standing-wave modes k, w(k) 

b=¢,, COs (7042/0) cos (7Ngy/b) COS (732/c) Exp (tt). Pet) 
The mode numbers 71, %9, %3 are positive integers, starting with unity, 
and ending where the wave number / 


ieee Ore AE 
an( 4 E+) are 


reaches the critical or ‘ cut-off’ value k,. 

In principle we could use this representation for any grain size, but it 
will be convenient to use it only if b, c are smaller than the lateral coherence 
widths. If the grains are larger this is not convenient because it is a 
redundant representation. It is evident that the standing character of 
the waves will not be noticeable if the grains very much exceed the 
coherence dimensions, but it can be shown that this will be the case even 
if they exceed them only a few times. Let us now consider 6, ¢ as the 
coherence widths, and B, C as the grain dimensions. First, it is clear that 
as soon as B>2b, C>2c there will be no restriction inside b, ¢ on the 
instantaneous configurations of the electric field. Yet, the grain walls 
could still be noticed by the waves sent out from b, ¢ which return into the 
region after reflection. But these reflections can be disregarded if the 
wave returns into b, c only after the interaction time. 

It will be shown later, and it is also evident, that in a grain smaller than 
the coherence region a single collision process will appreciably excite two 
modes in 21, 2g, N3 Simultaneously. Ina grain B, C, it will excite 2B/b, 2C/c 
modes in 7, 23, because these modes cannot be distinguished inside }, c. 
These modes then form a wavegroup, of approximate lengths b, c which will 
move towards the grain boundaries with the group velocity V,=dw/dk. 
With the assumption mentioned above this must be a small velocity, 
because of the strong dispersion. As an illustration let us use the law of 
Pines and Bohm for a free electron gas (degenerate). 


w=? +urk®, u®=8By/m + 5 ee 
where /, is the Fermi energy. This gives 
dw 6Hypk  6hkBy ._ Eyl 
dk Bmw SmB, By X 
Here EH, is the characteristic loss energy, and A=2z/k is the wavelength. 
E',/E;, is for most metals about 2/3, and the smallest (critical) ) is (in 
aluminium) about 104. Hence the group velocity is at most about 
6107 em sec”. 7 on the other hand, for a tungsten cathode, is about 
10-1 sec, hence the wavegroup will move at most about 60 Angstroms 
during the interaction time. This means that if the grain size is twice the 
coherence width, the standing-wave restriction will be absolutely un- 
noticeable inside a border 304 from the grain boundaries. In brief, as 


Vs 
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soon as the grains appreciably exceed the coherence region, we can 
substitute for the potential field inside b, c, instead of eqn. (2) the travelling- 
wave pattern 


TNX [Ng 
p=¢;, cos — exp {2x (2 “+ “= Vbexp (twt) tee (8) 


a 
in which the mode numbers 7», 3 run over positive and negative values, 
not excluding zero. It will be shown later that, other things being equal, 
the pattern 5 produces about four times the collision probability than the 
standing-wave pattern 2. 


§ 4. INTERACTION WITH VIBRATIONAL MopEs, 
AS A PERTURBATION PROBLEM 


Assuming first that the grain a,b,c is inside the coherent region, we 
describe the incident electron as a plane monochromatic wave, of finite 
duration 

W(x, y, 2, t)=const. exp {t(«gr—(E,/h)t} Ota 7 
where k)=27/X, is the wave number and H, the mean energy of the beam 
electrons. In the interaction this breaks up into a mixture of states 


V (dy, 2, =y,,(x, y, 2) exp {—1(H,,/h)t} 
with Ym= (abe)? exp {1(e+ Koy +k 32) }- 
If one subjects the wave numbers xy, «2, x, to the conditions 
Gk —27iM,, OK,=27Mz, CKs—27M, 
with integer m,, ms, m;—s, the %,, form an orthonormal set in a, b, c. 
One can then apply a well-known result of first-order perturbation theory. 


If at the time ¢=0 the electron is in a state O, the probability amplitude c,,, 
of the state m will grow during the interaction time according to 


dc .@ fo (opal is 
ea ms a) da dy dz 
dt izexp{ ( h, ye eae wim dbo uv oY 


where ¢ is the electrostatic potential in a, b,c. Assuming that the total 
of the transition probabilities is not large, the probability of the electron 
being found in state m after the interaction time 7 is 


e \2| 77 r@ re ce ; 2 
catia (7) | | | exp(—iQt)bn*dbby dt dx dydz|. . (6) 
h OO 0 


Here we have introduced the ‘ loss frequency ” 
Q=(Hy—#,,)h=LHyz|h 
and we will continue to measure energies in terms of frequencies. 

It is convenient and justified to consider the m-states as continuous, 
as the electron wavelength is always very small compared with the 
dimensions a, b, c. We specify the electron after the collision by its loss 
frequency , its polar angle @ and its azimuth ¢. Let the probability 
of the electron being found in the range dQ, dé, dé be 


p(Q, 6, ¢) sin 0d.Qd0db=c,*c, dn  . . . + (7) 
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p is the fraction of electrons per unit d92 and unit solid angle. dn is the 
number of m-states in this range. This is found to be 


abe 
sin 6x? dx dé dd 
(27)? 
abe m 
= —.— «sin dd dé dd. ey 
(Qr)e h* ° ( 
In the second equation we have gone over to polar coordinates in the 
x-space. In the third we have made use of the non-relativistic energy 


dn= se 3 ky dk, dks= 3 


(27) 


equation n= 2mE/h* 
and of Q=(Ey—E p)/h=h/2mM(Kg?—K pn”) 
hence d Q—= —(h/m)« de. 


In the last eqn. (8) we can replace « by xy without appreciable error. 
Combining eqns. (6) to (8) we obtain 


b 
p(Q, 8, d)= = = Koemn*Cm 
Sb m [ e\? 
—— i = a8 z = > 9 
One h GJ lhe re exp (—72t)if,,, th aes doa (9) 


This is a general formula for the ae and angular distribution of an 
electron beam which has interacted during a time 7+ with a potential 
field 6. We now assume that ¢ consists of a mixture of modes 

d=¢,, COS (70 42%/A) COS (TN yYy/b) COS (7N32/c) EXP (twt) tes 
with random phases, which are contained in the amplitude factors ¢,,. 
That is to say, in the mean 


Py $n*bm= 
if n and m stand for different sets of age MN, 3). Hence, carrying out the 
averaging process in eqn. (9), all cross-products wili vanish, and the 
result is a sum over the single modes. A longish but elementary calculation 


gives TACO Dare 2 
Ee | exp (i 2th, *dyjy dt dee dy dz] > 
“0.200100 
=7PLL L$, *$,1e][24][ro][rs] ee eS 
Ny Nz Ng 
where the square brackets are abbreviations for the following expressions 
: sin 7 iiss w)r \2 
sin ie —K,)a—7n,]|\2 = 2 
[n =( : ls ee ui 2) (Ko Ky)a , (12) 
2L(ky—K)a—77 | (Kg —Ky)a-+7N, 
rae (S 3(kyb—7N) \? Kb 2 
BN Eegb—705) Kob+an, } ’ a ee ae 


and similarly for [ns]. The first of these expressions is the ‘ energy 
resonance bracket’, the others are ‘momentum resonance brackets’. 
The first factor in (12) and (13) is the resonance bracket for a travelling 
wave, the second arises from a wave travelling in the opposite direction, 
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which together with the first produces a standing wave. It can be seen 
that when the fesy factor attains its maximum value, which is unity, the 
second factor is }. This arises from the fact that the effective part of the 
amplitude, corresponding to the first travelling wave, is one-half of the 
total, and the effect on the transition probability is one-quarter. 

In order to complete eqn. (10) we now only require the mean values 
¢,,*¢,. We obtain these from the consideration that in an electric wave, 
with wave number k, the electric field amplitude is & times the potential 
amplitude. The energy in the mode n is, in the mean, equal to the zero 
point energy $/w, and this again is equal to the electrostatic energy 
stored at the instant of maximum field amplitude. Thus 


— 125,76, = Hieo(h) eee eet: CLA) 


where k? is given by eqn. (3). Substituting this value into eqn. (10) w 
obtain finally 1 


28 b= samp eor*Z zelollnallnelits] + (18) 


where the summation is to be carried out over all the modes. from 
N4=N.=N,=1 up to the cut-off. 


§ 5. EVALUATION OF THE INTRAECTION PROBABILITY 


Figure | illustrates in physical terms the mathematical problem involved 
in the summation (15). For any electron emerging with wave numbers 
K4, Kg, K, the maximum of the contribution will come from those terms of 
the sum-over-modes which are nearest to resonance, that is to say for 
which the conservation of energy and momentum is most nearly satisfied. 
(A standing wave carries no momentum as a whole, but one of its 
components, which was effective in the collision, does.) 


Fig. 1 


tice point corresponds to a 

tum diagram in the X Y-plane. Each. lat Pp 

Meee code The hyperbolic regression line between the longitudinal and the 
transversal momentum arises from the combination of the conservation 


equations with the dispersion relation. 
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The momentum factors [7,], [72] can be imagined as resonance functions 
centring on the lattice points 71, m5. In fact these resonance functions 
overlap sufficiently for making the discontinuous mode-structure all but 
unnoticeable. In addition to the momentum factors there is, however, 
also an energy resonance factor [w], which has a maximum value where 
Q=w(k). It has been pointed out independently by Mott,* by Watanabe 
and Kanazawa,t and by Ferrell; that if the energy and momentum 
conservation conditions were exactly satisfied this would result in a 
unique deflection angle 6 for every loss Q. This can be shown as follows ; 

Conservation of momentum gives 


k2=ky?+K2—2xgk COS O=(Kky—K)?+ 4k qk Sin? 30. | i ae 
Conservation of energy gives, using the non-relativistic law 
Ko —Kk & (kg?—K?)/ 2x y= Q/v. ity ees 


Adding to these a dispersion relation, of the form 
w= w+ uk? 
putting Q=w and combining the above equations we obtain 


(caer B- (FS) =(S-a)@-Se. - 8) 


Uv 


where we have written «, for the transversal wave number of the electron. 
This is a hyperbolic relation between the loss 2 and the deflection angle @, 
which has been recently brilliantly confirmed by Watanabe,} giving, 
for Be, Al, Mg and Ge values of wu? in remarkably good agreement with 
the value u?=$H,/m derived by Pines (1954) for a free, degenerate 
electron gas. 
In terms of wave numbers, using eqn. (17), this gives another hyperbola 


2 Wrz bis 
P= {1} (ge 8 te 


which is shown in fig. 1. If the bombarding beam had an energy spread 
this hyperbola becomes a ‘ regression line ’, indicating the position of 
the mean energies. The energy distribution is, apart from the shift, the 
same at every angle @ and equal to the primary distribution. 

These somewhat rough considerations, based on rigid conservation laws 
which need not be exactly satisfied in thin films, give some useful hints 
for simplifying the rather complicated sum in eqn. (15). Using eqn. (17) 
we write this sum in the form 

w (sin 3{2—a(k)}7\2 (sin $(Qr7,—7N4)\?2 Qr, \2 
PE (So=amy) (Fae) (acta) Male - ey 
where we have introduced the electron transit time T,=a/v. 


It can be seen immediately that the energy resonance (the first factor), 


; —= 
* Private communication. 


} This information is taken from the article, still unpublished at the time of 
writing by D. Pines on ‘Electron Interaction in Metals’, to appear in 
Advances in Solid State Physics, Vol. I. Academic Press, New York, 1955. 
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is much sharper than the z-momentum resonance (the second factor), 
because in all thin films as used in experiments the wave transit time + 
is many times larger than the electron transit time Tx 

On the other hand it can be shown that in terms of mode numbers a 
the energy factor is much slower than the momentum factor, except for 
extremely thin films. If the mode number n, changes by unity, the 
argument of the sine in the second factor changes by 47, while in the 
first factor it changes by 

5t@o0dk  , Vark, 

°dk dn 8 ah 
We substitute here, using the free-gas dispersion relation 

V —»-=$(Ep/E,)(h/m)k, t=Th/kT and k,=Q/v=LH,/hw 
which gives 
dw dk 6 eh Vy _. Vp 
where Vy is the Fermi energy in volts. As an example, if Vp—10v, 
B=0-197 and T=2500°c, this gives 
$r(1-72 x 10-7/a). 

That is to say the energy resonance becomes as sharp as the momentum 
resonance in terms of the mode number n, only at thicknesses a of 
17 angstr6m. Thus for thicknesses of more than 30-604 we are justified 
in considering |w]| as a slow factor, and to put it before the summation sign 
in eqn. (15). This means that for all except very thin films the energy 
distribution will be the same as the primary distribution, except for a 
shift as calculated from eqn. (18), which is a function of the angle 0. 
(The sum will still contain 2, but will be extremely insensitive compared 
with the factor before the sum.) 

For very thin films a more accurate calculation would be needed, but 
this is hardly worth while because it would be extremely difficult to check 
experimentally the slight departures of energy distribution in function of 
the angle @ which it entails. Moreover, these small effects vanish altogether 
in the integrated probability, because [w] is so much sharper than the 
rest in terms of 2, especially for thin films. 

Putting [w] before the summation sign, we are left with 


pel TE sin $(w97,—74) \? Oye, 

>| o/{*(i é B u: =)} | ( 2(@7-—771) ee allel: 
| ree Tea 
Here we have replaced w in the first factor and 2 in the next two by @p, 
as this makes no significant difference. 

Consider now the angular distribution. This is entirely contained. in 
the factors [7] and [ng] as only these contain x, and x3. We can however 
simplify matters by replacing the first factor under the sum by 

Wo 


(71/0)? + ko?+ Kk” 
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that is to say replacing the transverse mode numbers 7%, 13 by ee 

resonance values. This again is a slow factor, and can be put ere the 
7 7 

summation signs of the sums in n, and n3. We are now left wit 


“el SIN }(WoT,—7N 1) Wo Te q) a a(n 
Cewseeneordt ea, St ot ae 


which is the product of three sums. The last two can be written in the 
form 


2 [N9]=F (kb/7), 2 [N3]=F (Kgc/77) . 3 Ses 

in terms of the function ‘ 
N /sin }n(x—n)\? { x ) (94 
f0)= 2 (ey) ean) 


This function is illustrated in fig. 2. It starts at zero for ~=0, rises rapidly 
near x=1 and attains a maximum of 0-5 at ~=2. It then approaches, 
with insignificant fluctuations, the asymptotic value 0-481. 


Fig. 2 
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The function F(x) and its components. 


It might appear from this that the polar diagram falls to zero at zero 
angle 6, moreover that there is a cross-shaped gap in it at é=0 and d=4r, 
parallel to the axes. But these are just covered up by the diffraction 
which we have neglected when representing the emerging electron by a 
plane wave, instead of a wave-bundle which has zero value outside bres 
It is therefore more accurate to substitute in the eqn. (22) the asymptotic 
value 0-481 of F(a), up to the cut-off, and it is seen that the angular and 
energy distribution is given by a sum of terms of the form 

[o] - 
(rma) (eae, ae) 
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For films which are not too thin one can moreover replace (7n,/a) by the 
constant value 2/v, which follows from the X-momentum resonance 
(cf. Ferrell 1955). This is illustrated in fig. 3. It appears that this type 
of angular distribution is in good agreement with experiments. 

We now calculate the integrated, total probability of excitation, P. 
The energy dependence is contained in the [«] factor, and this is so sharp 
that we can integrate between infinite limits, using the formula 


fo bide [" (Se aot. Pmen(26) 


Integration of the angular factor 25 over all directions, up to the cut-off 
gives 


(Ke) 6 dé 27 (* kK, dk, 7 k, 
|, (7n,/a)P?+KP = (7n,/a)P+K?P Ke sloe( =a) oe 


0 


Fig. 3 
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Angular distribution of the current density i per unit solid angle. 
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Here we have integrated up to a maximum wave number /, starting from 
zero transverse wave number, in accordance with the remark made above 
that diffraction makes the drop-out of the modes n,=0 and n,;=0 
unnoticeable. 

Collecting all factors, from eqn. (15) on, we have now 


P= ||| p(2,0,4)a:2040dd— eee (<) wor F) 


; / 2 p 
(ee Z(WoTe~-71) oT! log ies 
x 
1 $(woTe—TN 1) WoTe TOM ee 


; (28) 


where we have substituted «.=mv/f. In the last remaining sum we now 
introduce the new variable 


L= WT ,|[T= Wp t/10=4/4,,. Ae me) 
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This is a very convenient measure of the thickness, in units of the 
thickness a, in which the electron ‘ sees ’ just half a cycle of wy. This is 


A, =| wy=nBch/eV,=6-12X10-*B/Vzem. ee (30) 


where V, is the characteristic loss in volts. For 10 kev electrons, S=0-197 
and aluminium, V,—14-7 v, this gives a,=80A. 
We now introduce the function 
NV sin d(a—n) a \? 
(x)= 2X n {| ———_ — 31 
G(x) log n( ery ary: (31) 


iT 

which is illustrated in fig. 4. In terms of this function and the previously 
defined F(x) (eqn. (24)), we can now write the last sum in eqn. (28) in 
the form 


ka Cpe. 
Sing} ton( = ( i a 


The functions F(x) and log «F(«)—G(z). 


The function log xF(x)—G(x) is also illustrated in fig. 4. It goes asymp- 
totically to zero after one negative and one positive excursion. 
Our final result is, therefore, for coherently illuminated grains 


__ (0-481)? 1 /c Pera oe 
P= 7h isi( $) eno (, ae Beer BEE 


a 


where we have used fc/e?=137. Instead of the cut-off wave number k 
we have used the corresponding critical wavelength \,=2z7/k The 
function Q is illustrated in fig. 5. : te. 

It is seen that the cross section or collision probability for a collective 
interaction behaves indeed very differently from the corresponding 
quantity in an atomic process. Instead of being linear in function of the 
thickness a, it starts first very flat, and then approaches rapidly a constant 
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asymptotic value which it overshoots slightly at 2a,. This asymptotic 
value is 

(0-481)3 1 /e _ WoT 
P.= (; (o97) log (a,/42,)=6-5% 10-9 27 log (ag/4A.). (34) 


where we must substitute 


ayrmeV kl =r 1LGOOV IT. 2 1 5. . (35) 


As an example take aluminium. The maximum deflection angle 0, 
which can be roughly identified with the cut-off, observed by Watanabe 
with 25 kev electrons was 1-2 10-%. This corresponds to X,=9A4. With 


V,=—14-7 v and a tungsten cathode wyr=214. This gives for 
V=10kv, B=0-197, a,= 804, P=0-204, 


V=25kv, B=0-304, a,=1264, P=0-148. 


Fig. 5 
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The collision probability function Q, in function of the thickness and of the 
cut-off wavelength A,. 


Before discussing this and other consequences of eqn. (33), we calculate 
the corresponding quantities for incoherent illumination. 


$6. INTERACTION WITH GRAINS SIGNIFICANTLY LARGER 
THAN THE COHERENCE WIDTH 
In this case we have essentially incoherent illumination. As explained 
in § 3, each elementary excitation process produces travelling waves of 
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the type as in eqn. (5). Nothing changes in the previous calculation up 
to eqn. (13), which we must now replace by 


ara sin 4(kyb—27N9)\? 13’) 
pral'= ( epaaae) ee 


and correspondingly eqn. (24) must be replaced by 
N (sin 4n(a—2n)\? re. 
ita —_———_]} . “a. =, ee 
AG >i In(@—2n) ae 
But this function is exactly unity, independently of x, if V goes to infinity, 
and if N is finite it is still appreciably unity for —N<x<WN and zero 
outside. Hence in the incoherent case nothing is changed in the energy 
or angular distribution, but the previous factors 0-481 in eqn. (33) must 
now be replaced by unity. This means that in the case of incoherent 
illumination the probabilities P and P., in eqns. (33) and (34) must be 
multiplied by 1/0-4812?=4-32. 


§ 7. COMPARISON WITH THE EXPERIMENTS OF GABOR AND JULL 


In these experiments, which have been only partly published so far 
(Gabor and Jull 1955), ‘ stratified ’ layers of metals were built up of thin 
layers, with some foreign material between. In one experiment for 
instance six layers of gold of 25-304 each were evaporated on top of one 
another, with thin layers of perhaps 10-204 of ‘ carbonaceous con- 
tamination ’ between, which are easily produced in any continuously 
evacuated system by electron bombardment, which cracks the oil and 
grease molecules which are always present in the residual gas. It was 
found that the 16 ev gold line was just perceptible in the layer with 
six strata equivalent to 170A of solid gold ; about an order of magnitude 
weaker than in a solid gold layer of 110A thickness. 

These experiments undoubtedly indicate that the initial rise of P is 
stronger than linear, the difficulty is only that the evidence is too strong. 
For 10 kev electrons, as used by Gabor and Jull, and V;—16 vy, a, is 
about 75A hence the six layers correspond to about z=0-38 each, mene 
the single gold layer of 110A corresponds to x=1-5. As seen in fig. 5 
such a strong effect could be expected only if a,/}A, were about 2-5, Soi 
corresponds to a cut-off wavelength of about 604. The first interpretation 
put forward by Gabor and Jull was, accordingly, that the critical wave- 
length must be much longer than expected on the basis of the theory of 
Pines and Bohm. In view of the experiments of Watanabe, which indicate 
cut-off wavelengths of rather less than 104 we must now explore 
alternative explanations. ; 

One possibility suggested by the foregoing theory is that perhaps the 
grain size in the 25-304 gold films was much smaller than that in the 
foil of 1104 so that the first were coherently illuminated and the other 
incoherently. Ifthis were the case, one could just reconcile the experiments 
with a cut-off wavelength of about 104. This would Boprokimea 
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correspond to a curve for a,/$\,=15 in fig. 5. This gives Y=0-10 for 
+—=0-38 and Q=1-27 for x=1-5. Hence six layers of x—0-38 would give 
an effect of 0-61, but if the solid foil was incoherently illuminated this 
would have to be compared with 4-32 1-27=5-5, which is nine times 
stronger. But this explanation does not agree with direct observations 
of the grain size in electron micrographs and in diffraction diagrams. 
The electron micrographs showed indeed that the granules of the stratified 
layer were smaller, but they were still of the order 10004 which is much 
more than the coherence width used in these experiments. Electron 
diffraction diagrams even indicated that the stratified layers were better 
crystallized, which is quite likely in view of their more prolonged 
bombardment. It appears therefore that this interpretation must be 
dropped. 

There remains only one likely explanation, and this is that it is too 
naive to think of layers of the order 25-30 A thickness as if they had the 
properties of matter in bulk, as regards collective oscillations. If one 
imagines that there is on the surface of the metal foils a ‘ disturbed ’ 
surface layer of only 3-4 atomic layers thickness, so that the standing waves 
have these as boundary, this would mean that we have to consider the 
effective thickness of these layers as perhaps 15-204 less than the nominal 
thickness and the experiments are readily explained. 


__ § 8. ConcLustons. NEw Cructau TEsts 

Consider again eqn. (33) for the probability of exciting the first-order 

characteristic loss 
: -4 a ; 

Pm He at (5) (WoT) Q (= ; =) cee illumination. . (36) 

This function has so many features in which it differs from the probability 
of exciting single-electron transitions, that it ought to be possible to decide 
about the ‘ collective’ or ‘ single-electron ’ nature of any characteristic 
loss beyond a doubt. We have already discussed the non-linear start for 
small thicknesses. Another striking feature is the asymptotic value of P, 
which is reached at about a= 2a,, that is to say at a thickness of the order 
100-1504 under usual experimental conditions. This is easily understood. 
Once the thickness is of this order, any single collision act will excite 
significantly 2-3 modes in the X-direction. This number remains constant 
if the thickness is increased ; the momentum resonance becomes sharper 
instead of stronger. The reason is that the interaction time is determined 
by the train length instead of by the foil thickness. 

We have seen that P., for coherent illumination is of the order 0-15-0-2, 
while for incoherent irradiation it is of the order unity. In this case of 
course first-order perturbation theory is no longer accurate. In reality 
P would increase, once it has reached appreciable values, slower than 
the Q-function, because the primary electrons are used up, nor would it 
go beyond about 0-5 because 2d, 3d, ... order transitions would become 
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important. From the fact that some experimenters have observed 
simultaneously lines of five orders in Al, we must infer that in these 
experiments the illumination was incoherent. 

The eqn. (36) suggests two crucial tests. The first arises from the 
factor wy7 which is inversely proportional to the primary energy-spread, 
i.e. to the temperature of the cathode. With an oxide cathode, 7’=1000°, 
P should be 2-5 times stronger than with a tungsten cathode, T7=2500°. 
This means that, using the same foil, weaker lines, if they are due to 
collective interaction, should be increased in intensity 2-5 times, while in 
the case of stronger lines the higher orders should increase very noticeably. 
It appears that all experimenters have used tungsten cathodes, as is usual 
in demountable vacuum devices. Ifthe hypothesis of collective interaction 
is right, it ought to be possible, using oxide cathodes or some other methods 
of producing electron beams of more uniform energy, to produce the 
collective part of the loss-spectrum almost pure, much stronger than 
hitherto observed, and to separate it from loss-lines of other origin. 

A second crucial test is suggested by the large difference between P for 
coherent and for incoherent illumination. Most experimenters have 
probably operated with rather large grains (well-crystallized foils), and 
not too well collimated electron beams. With grains of the order of a 
few hundred angstroms it ought to be possible to notice a strong diminution 
of the characteristic losses of collective origin if one passes from, say, 
critical illumination to a well-collimated beam, with divergences of less 
than 10~¢ radian. . 

It may be noted that this theory and the proposed crucial tests apply 
only to sharp lines. This is the consequence of our assumption that the 
collective oscillations are undamped, up to the critical wave number, 
where they are suddenly cut off. If they are appreciably damped, as is 
probably the case in insulators, and also in most heavy metals, the theory 
must be revised, introducing complex values of w, to account for the 
‘broad lines observed in many experiments. 
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APPEND DX 
THE EQUIVALENT TRAIN LENGTH OF ELECTRONS 
Emirrep By A Hot CatHopE 


:; Itz, is the initial velocity of the electrons in the x-direction, the current 
in the limits v9, dvy is, at the cathode, proportional to 


EXp (— {Mv _7/kT)v9 dvg. 
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Adding an energy EZ, corresponding to a voltage drop E 
istribution in ; e, th ? 
distribution in the beam current becomes . P Ho peer ey, 


exp |—(H—E,)/kT] dE. 
The mean energy of the electrons in the beam is £y,-kT'. This distribution 
is shown in fig. 6. 

Compare this with the energy distribution corresponding to a chopped 
wave train of frequency w, and duration 7. To obtain this we need only 
calculate the absolute square of its Fourier transform, going over to 
energies by H=hw. This gives 
(= eee _ (sin ((E—E£y—kT)t/h}\? 

4(w—wy)t = 3(L—Hy—kT)r/h ) 
where we have identified wh with Hj+kT. Figure 6 shows this distribution 
for three values of 7, $rf/kT, wh/kT and 2ch/kT. 


is Q / 2 —+(E-E\kT3 


Energy distribution in an electron beam emitted by a hot cathode, compared 
with three chopped monokinetic pulses of different durations. All 
areas are equal. 


Tn order to make the best choice we make the comparison in yet another 
way, by comparing the chopped monochromatic train with the ‘ time- 
representation ’ of an electron beam with Maxwellian probability distri- 
bution. We obtain this formally by taking the square root of the 
Maxwellian density in £, and then taking the absolute square of its 
Fourier transform, which is of the form 

1 
1+ (2kTt/h)? ~ 


The physical meaning of this formal representation is that the electron 
beam issuing from a hot cathode of temperature 7' could not be distinguished 
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from a source which emits probability pulses of this shape. This pulse shape 
is compared with chopped pulses of durations r=$7h/kT, 7wh/kT and 
2rh/kT in fig. 7. 


Fig. 7 
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The probability pulse corresponding to a Maxwellian beam, in time represen- 
tation, compared with three chopped waves of different durations. All 
areas are equal. 


Testing the fit by the rule of least squares in both representations, the 
conclusion is that r=7h/kT' is the best compromise. The energy repre- 
sentation appears to favour slightly larger values, while the time 
representation favours smaller ones. The compromise solution t=7h/kT 
has been adopted in the text. 
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Il. The Optical Model Applied to the Elastic Scattering of Nucleons by 
Various Light Nuclei 
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H. H. Wills Physical Laboratory, University of Bristol § 


[Received September 14, 1955] 


SUMMARY 

The Optical Model of the nucleus has been applied to the elastic 
scattering of 9-5 Mev protons by He, C, N, O, F, Ne and A, and 14 mev 
neutrons by C, N and O. In general, fair agreement has been obtained 
for a square-well potential, but the results suggest that for higher energies 
and heavier nuclei a shaping of the well is necessary. By making the 
nuclear radius a variable parameter, instead of 1-45 A3x 10-33 em, 
better agreement could probably be obtained. The mean free path of 
nucleons of these energies in nuclear matter is found to be larger than the 
nuclear radius. No simple relation was found between the depth of the 
potential-well and the atomic weight. 


§ 1. INTRODUCTION 
THE Optical Model of the nucleus (see e.g. Feshbach e¢ a]. 1954) attempts 
to reconcile the independent-particle model with the strong-interaction 
model. At the same time, it allows the scattering of a proton or neutron 
by a nucleus to be treated as a two-body problem. The nucleus is 
represented by a potential-well of complex depth, and of radius R, 
such that in the simplest case of a square well 


V= — ( Voti Wo); for r<R. 
The Coulomb effect for protons is represented by 
V=Zetir, for ¢>f. 


It has been shown by Chase and Rohrlich (1954) that the square-well 
potential cannot give a good fit to the angular distribution of 18 Mev 
protons elastically scattered from aluminium. They conclude that 
agreement may be better for nuclei with smaller Z. However, the work 
of Britten (1952) had already demonstrated that this simple theory does 
not fit for beryllium bombarded by 31:5 Mev protons. Tt appeared 
therefore that an application of the theory to the scattering of lower 
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energy protons (9:5 Mey) by the lighter nuclei would be of interest. At 
this energy, the wavelength of the incident nucleons is of the same order 
as the radius of the nucleus and therefore the effect of the edge is not so 
serious as it is for higher energies. Also the calculations are somewhat 
easier, especially in the absence of an electronic calculating machine, 
since only a few values of L (usually L=0-> 4) are significant. 

For reasons given by Adair (1954) from low energy neutron data, the 
expected value of Vy is in the region of 42 Mey. Feshbach et al. (1954) 
give theoretical curves of differential cross sections for | Mev neutrons 
elastically scattered from elements of A=40 to A=260. They use 
V=—42(1-+70-03), and the agreement is, in general, remarkably good. 

A two-body interaction model leads to a more accurate representation 
of the scattering of nucleons by even nuclei than by odd nuclei. With 
the exception of nitrogen and fluorine, the elements studied here are 
even-even. It is of considerable interest, therefore, to compare the 
results of the application of the Optical Model to these two nuclei (N, F) 
with those obtained for the other nuclei. We can expect a small value of 
W, for even nuclei, which have few low lying levels. 


§ 2, APPLICATION TO OXYGEN 


Oxygen is a closed shell even—even nucleus and it was for this element 
that the properties mentioned above were demonstrated experimentally. 
Thus it has been shown (Ajzenberg and Lauritsen 1952) that the value of 
the reduced width indicates the validity of the two-body model for low 
energies. The value of the resonant energy (H,) for neutrons can be 
related to the radius R and potential V9 of a square well by the equation, 
E,=—(Vo+E#,) cot? (k’R), where k’=V2M(V)+#,) since only waves of 
L=0 are involved. This gives, for R=3:65x10-%cm, the value 
Vo=42 Mev. At higher energies, however, other values of L are involved. 
Hence the criteria for fit allow, in principle, a less ambiguous assignment of 
V and Rk. On the other hand, this is complicated by the introduction of 
Wo, which in the case of oxygen is conveniently small. In an earlier 
communication (Fujimoto and Hossain 1955) some of the results for 
oxygen were reported. This element was chosen for the reasons given 
above for an extended study of the variation of the differential cross 
section curve with the three parameters; R the radius of the nucleus 
and V» and Wy, which determine the complex potential V=—(V,+iW, ) 

The expression given by Blatt and Weisskopf (p. 336) can be aa 
rearranged using the same notation to give the more convenient form 


do(d) _ [ Ze* ed [boss 2 
AQ = | Tyr c08ee (¢/2)— i ema) cos wy, sin 6,Pz (cos 6) | 


] ce } 2 
= a) 2 CL+ 1) sin %, sin 6, P,; (cos 6) | : 


where p= 2(o,—09) +874 2y log, (sin $/2). 
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The value of R for the present calculations was given by 
R=r, A3x10-% em, where A—Mass Number. Following Feshbach 
et al. (1954) the value of ry was taken as 1-45 giving R=3-65 x 10-33 em 
for oxygen. The contributions from terms involving L>4 were shown 
to be negligible in the case of oxygen. It was found that the shape of 
the curve depended critically on the value of V, as shown in fig. 1. As 
V, was varied from 32 to 42 Mev with W,=0 the shape showed a marked 
change at V)=38 Mev. This was occasioned by the magnitude and sign 
of the f-wave (L=3) component near 6=70° which first cancelled out the 
other terms for V)~36, and then at Vj;=38 swamped the L=0, 1. 2 
contributions and produced a subsidiary maximum. 


Fig. 1 
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Theoretical proton scattering curves. o(f) is in mB/steradian and ¢ is the 
c.m.s. angle (as for all figures). 
This rather complicated behaviour is also presented by ae elements 
i Ave i itrogen). owever, in 
for certain values of R, V, and W, (cf. fig. 5, for nitrogen) ae 
; tage: 
attempts to fit the experimental results, this leads to no real difficulty. 
Preliminary values, for neon in particular, sometimes gave curves near 


an anomalous region. 
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" ‘The introduction of W, tends to flatten the curves, in the sense that it 
increases o(¢) at the first minimum and reduces the first maximum. 
With the exception of near-critical values of Vg it does not alter very much 
the positions of the first minimum and maximum. This was found to be 
true for all the elements studied, and the curves for nitrogen (fig. 5) 
show both the typical and the anomalous dependence on Wo. 

Two criteria were found to be useful for determining the ‘best’ fit. 
The angular position of the first interference minimum was adjusted 
by suitable variation of V,. In most cases this also served to locate the 
first maximum at approximately the correct value of ¢. The other 
criterion was the value of o(f) at the first maximum and this could 
usually be obtained by adjusting W5, remembering that increasing Wo 
was always found to lower o(¢) at this maximum. It was expected that 
the departure from the experimental values would be greatest for the 
higher values of ¢. This justified the relative neglect of the position and 
absolute value of the second minimum. 

The theoretical curves giving the closest fit for oxygen are shown in 
fig. 2; the experimental points are from Burcham ef al. (1953), with 
certain additions communicated privately. The curve for )=—(42-+70) 
is also shown because of the importance of this value of V, in the interpre- 
tation of neutron data at 1 Mev. There is a series of values of V) from 
30 to 36 Mev, with a suitable value of W, from 1 to 2 Mev, giving curves 
in very fair agreement with experiment. The values near V,—36 Mev 
give a slightly better value of dé for the first minimum. The absolute 
value of o(¢) in this region could most probably be improved by tedious 
adjustments of V, and W,, since V, is becoming very critical. As already 
indicated, at Vj)=36 Mev the maximum contribution to the scattering 
comes from proton waves of angular momentum Z=3. This suggests 
that the states of F near 10 Mev are predominantly f-states and hence 
have spins of 5/2 or 7/2, thereby showing that the one-particle model is 
in this case a sufficiently good approximation for explaining the energy 
states concerned. 

The application of the Optical Model to the elastic scattering of 14-1 Mev 
neutrons by oxygen, as already reported (Fujimoto and Hossain 1955) 
showed that the values V=—(32+i7) and V=—(36+i5) give a fairly 
good fit to both the differential and total cross-sections. However, the 
experimental results (Conner 1953) for the differential cross sections are, 
as yet, insufficiently precise to distinguish between the possible values of V. 
The values of total cross section are of greater interest. As shown in 
table 1, the value V==—(32+78) gives the better ratio of the absorption 
and elastic cross sections for a total cross section equal to that obtained 
experimentally (1-6 barns). Nevertheless, the ratio (0-70) is not in good 
agreement with the experimental value (1-28). The best ratio for cae 
42 Mev gives a poor value of total cross section, but still a tolerably good 
fit to the angular distribution. This value, V,—42 mev, is of further 
interest in that it gives a bound 2-s state for 170 at a binding energy, for 
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Fig. 2 
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Experimental points and theoretical curves for 9-5 Mev protons elastically 
scattered by oxygen. 


Table 1. Cross sections for the Scattering of 14 mev Neutrons by 
Oxygen (in barns) 


Square- Well Elastic (nett 

: . ption Total 

ee Salm cross section eross section 
30+725 1-21 0-59 1-80 
32472 1-68 0-43 2-11 
32475 1-16 0-61 1-77 
32476 1-08 0-63 1:71 
32+78 0-93 0-66 1-59 
36475 1-01 0-54 1-55 
36+76 0-96 0-58 1-54 
40+25 0-70 0-44 1-14 
42425 0-73 0-45 1:17 
42418 0:66 0-59 1-25 

Experimental 0-7 0-9 1-6 
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the well known 0-875 Mev level, of 3-23 Mev in exact agreement with the 


experimental value. 
The effect of varying V, or W, can be seen from the tables. General 


rules can easily be drawn up, but certain exceptions have been found. 


§ 3. CARBON 
Using the same criteria of fit for protons elastically scattered from 
carbon, it was found that a square-well potential with V)~44 Mev and 
W,=4 Mev gave the best results (fig. 3). In this case V)—42 Mev is a - 
tolerably good fit, but the position of the first minimum suggests that 


Fig. 3 
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Experimental points and theoretical curves for 9-5 Mey protons elastically 
scattered by carbon. 


higher values of Vy may give a better fit. The curves are for L <3 
higher values of L making negligible contributions. The experimental 
points given here are those of Burcham e¢ al. (1953) with certain additions 
communicated privately. The results of Fischer (1954) have been com- 
ae and sieees e ee a somewhat smaller cross section. For ¢6>60° 
ney are we e V=— LO 7 ‘| 
ta laa oe (42+76°3), but for 6<60° the experimental 
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Experimental points and theoretical curves for 14 mev neutrons elastically 
scattered by carbon. 


Table 2. Cross sections for the Scattering of 14 Mev Neutrons by Carbon 
(in barns) 


Square- Well Elastic Absorption Total 
Potential scattering cross section cross section 
(Mev) cross section 
32+78 0-82 0-49 1-31 
37+16 0-87 0-48 1-35 
3717 ett. 0-54 1-3] 
42475 0-85 0-56 1-41 
42476 0-76 0-56 1-32 
42418 0-61 0-57 1-18 
44-75 0-81 0-54 1-35 
44-16 0-73 0-56 1-29 
47+-15 0-69 0-49 1-18 
Experimental 0-72 +0-04 0-60 0-03 1-32 40-02 | 
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The larger values of W, for this element, compared with that for oxygen, 
may well be explained by the presence of more low lying levels in carbon 
than in oxygen, providing a greater total inelastic cross section. The 
thresholds for both (p, n) and (p, pn) reactions are above 9-5 Mev, and 
the cross sections for (p, «) and (p, y) are probably very small. The values 
of ry for carbon, nitrogen and argon are slightly different from that for 
oxygen. This was occasioned by a small simplification of the calcula- 
tions, and is known to have a negligible effect on the theoretical curve. 

Preliminary calculations have been made in an attempt to fit the 
experimental results of Wright (1953) for the scattering of 30-6 Mev 
protons by carbon. The radius used was 3-320 x 10-1 cm and values of 
L <4 were considered. It was possible to give a fairly good fit up to 
6=90°, for V=—(42+76-3) but the positions of the higher angle minima 
were at much smaller angles than the experimental values. The cross 
sections for 6>140° rise steeply to values near 180° which are of the order 
of 100 times the experimental values. Both of these divergencies confirm 
the work of Chase and Rohrlich (1954) and indicate the greater need for 
shaping the well at higher energies, as done by Woods and Saxon (1954). 

The differential cross section for 14 Mev neutrons elastically scattered 
from carbon is available only for 6 <55° (Coon, private communication). 
This is insufficient for a determination of the best values of V, and Wy. 
However the further criterion of total elastic and absorption cross sections 
can be applied. Table 2 gives the values for a series of potentials, and 
the experimental results. Combinations of V, and W, are given which 
lead to a total cross section of about 1:32 barns, i.e. the experimental 
value. As Vy is increased from 32 to 47 Mev, W, must be decreased from 
8 to 4 Mev. At the same time the absorption cross section increases and 
the elastic scattering cross section decreases, each approaching the experi- 
mental value. It is seen that the most probable value of V is —(42+76). 
Figure 4 shows the corresponding differential cross sections, and the 
experimental points. 


§ 4, NITROGEN 


The experimental values of the differential cross section for 9-5 Mev 
protons elastically scattered from nitrogen (Freemantle et al. 1954) are 
shown in fig. 5. The points do not appear to lie on a smooth curve and 
the values rise steeply at backward angles (6>150°). It is possible that 
this is associated with the odd-odd system of nucleons, which can be 
expected to produce a complicated configuration of the nucleus. The 
Optical Model gives no more than a moderately good fit. The best curve 
is that for V=—(42+74-2). Figure 5 includes the curves for V o=42 
with W,)=0, 2-1, 4-2, and 6-3 mev showing the characteristic variation 
of the shape as W, is decreased from 6:3 to 2:1 Mev. The curve for 
W)=0 shows the anomalous subsidiary maximum near d= 60° occasioned 


by a large contribution for L=3 (ef. fig. 1 where this occurs atV,—38 
with W,=0 for oxygen). 
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The differential cross section for 14 mev neutrons elastically scattered. 
from nitrogen has been measured by Smith (1954) and the total elastic 
and absorption cross sections are given by Coon et al. (1952). Total 
cross sections are easier to calculate than differential cross sections and 
the results for certain potentials are given in table 3. No really good 
agreement could be obtained, the best values being for V=— (36478) 
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Experimental points for 9-5 mev protons elastically scattered by nitrogen 


The series of theoretical curves shows the variation with W, as discussed 


im the text. 


Difficulty was experienced in obtaining a sufficiently high value of the 
total cross section consistent with an approximately correct ratio of aie 
to absorption cross section. It was also difficult to fit the eee 
cross section curve, as can be seen in fig. 6. The general pattern o 

maxima and minima is of course reproduced, but the absolute values are 
far from correct. Thus it is not possible to make a sensible comparison 
between the best values of V for protons and neutrons. 


This also may be 
due to the odd—odd system of nucleons. 


On the other hand, the paucity 
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Experimental points and theoretical curves of 14 Mev neutrons elastically 
scattered by nitrogen. 


Table 3. Cross sections for the scattering of 14 Mev neutrons by Nitrogen 
(in barns) 


Square- Well Elastic ; 
Potential scattering Absorption Total : 
(Mev) Grote Bonan cross section cross section 


324-78 0-89 (5) 0-58 (5) 
36+i8 0-86 0-63 (5) 
39+-78 0-80 0-62 
42+-75 1:03 0-49 
424-18 0-71 0-56 
48 +15 0-95 (5) 0-30 
48 +78 0-57 0-49 
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mt 


of differential cross section data for neutrons of this energy does not allow 
the conclusion to be drawn that nitrogen is unusual for both neutrons 
and protons. It is possible that the square-well Optical Model does not 
give a good representation of the differential cross section curve for 
neutrons of this energy region—but in the light of its fair success for 
protons this would be a surprising result. 


Fig. 7 
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Experimental points and theoretical curves for 9:5 Mev protons elastically 
scattered by helium. 


§ 5. Hetium, Neon, FLUORINE AND ARGON 


Angular distributions for 9-5 Mev protons elastically scattered from 
these elements have been obtained experimentally. The values for 
fluorine have been kindly communicated, prior to publication, by 
Dr. D. J. A. Prowse. For the other elements the results used here ae cate 
of Freemantle et al. (1954) with certain additions communicated payate y- 

Other data are available for helium, and those obtained by Kae 
(1952) have been subjected to a detailed phase-shift analysis by Dodder os 
Gammel (1952). Nevertheless it was felt to be worthwhile to use the 
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Optical Model for this element also, in order to discover the limits within . 
which it could be applied. The results are shown in fig. 7. It was found 
that contributions arising from L >3 were negligible because of the small 
size of the nucleus. The fit is better than might have been expected 
especially considering the relative crudity of the model which neglects 
inter alia spin-orbit coupling. The value of Wy was taken to be zero, 
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Experimental points and theoretical curves for 9-5 Mev protons elastically 
scattered by neon and fluorine. Note that the o(¢) scale for neon is 
displaced upwards by half a cycle. | 


thereby assuming a zero absorption cross section. However, we tenta- 
tively suggest that the fit may be improved by introducing a value of 
W, to allow for the formation of the compound nucleus ®Li and then 
calculating the ‘ compound-elastic’ scattering. In this particular case 
the only channel open to the protons is that leading to the formation of an 
alpha particle. It is also possible that a larger radius could be invoked 
to produce a larger cross section and therefore a better fit. | | 
Neon and fluorine yield almost identical experimental values of angular 
distributions, and are shown for comparison in the same figure (fig. 8). 
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The o(¢) scales are displaced to show the agreement with the theoretical 
curves. This agreement is quite good, except near the second sana 
for potential values of V=—(34-5+76-8) for neon and ee es, 
for fluorine. The higher value of W, for neon, compared with nitrogen i 
consistent with the available energy levels and the sum of ie es 
sections, which is fairly high. In the case of fluorine however. the ee 
of W, might be expected to be larger than that found here It would 
appear from these results that the square-well Optical Model applies 
equally well to the even nucleus neon and the odd nucleus fluorine ra Ain 
interpretation of the scattering of 9-5 Mev protons. The close resemblance 
of the experimental curves for these two nuclei is very striking. 


Fig. 9 


Experimental points and theoretical curves for 9:5 Mev protons elastically 
scattered by argon. 

Calculations have also been made for 9-5 mev protons elastically 
scattered from argon. The results are shown in fig. 9. For the radius 
used here, 4:80 10-13 cm, the value of o(¢) is too high, suggesting that 
it might be more nearly fitted by using a smaller radius. However, it is 
possible to fit the angular position of the first minimum and maximum 


by V=—(42+ 14-2). 
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§ 6. Discussion 


The present results indicate that the scattering of low energy (9:5 Mev) 
protons by the lighter nuclei can be moderately well represented by the 
Optical Model of the nucleus, using a square-well potential as a first 
approximation. It is however confirmed that scattering of the same energy 
protons by somewhat heavier nuclei (e.g. argon) cannot be so well accounted 
for by such a simple model. This agrees with the conclusions of Chase 
and Rohrlich (1954). The work of Britten (1952) on 31-5 Mev protons 
scattered by beryllium appears to be confirmed by the preliminary caleula- 
tions for 30-6 Mev protons on carbon, showing that higher energies also 
require a more complicated model for a satisfactory account of experi- 
mental values. 

It is quite possible that better fits could be obtained by a systematic 
investigation of the influence of the radius R on the calculations. In the 
case of oxygen, preliminary calculations showed that a radius of 
4-4 10-13 em gave much too great a cross section. Also a change from 
3-65 to 3-7510-33 cm made very little difference to the shape and 
absolute cross section for a given potential. In this connection it is 
noticed that the formula R=(1-241/8-+-0-7) x 10-18 em recently proposed by 
Weisskopf (private communication) gives values of the radius which would 
lead to better fits, as discussed above. Thus for helium this formula gives 
R=2-61x10-% cm compared with the value 2-3 10-18 cm used here, 
and for argon R=4-8x10- cm instead of the value 4:95 <10-% cm. 
In the latter case the change in the differential cross section curve would 
be very small. 

For the elements considered here, the mean free path of the incident 
protons in nuclear matter lies between 4-410-!% cm (for neon) and 
8-2 x 10-18 cm (for carbon), except for oxygen where it is somewhat larger 
(15-1 10-13 em). These values are calculated from the relation 
4-5(L+- V)1/?/W, (see Chase and Rohrlich 1952) using the best values of 
V, and Wy. No value was found smaller than the nuclear radius. 

It is possible that a particular element is characterized by a certain 
value of Vy but this may vary with the energy of the incident nucleons. 
This points to the need for more experimental values of the differential 
cross sections of elastically scattered neutrons at various energies. The 
calculations for neutrons are much simpler, and the absence of Coulomb 
scattering allows comparisons of total cross sections. 

The work reported here indicates the limits of application of the square- 
well Optical Model. A closer approximation to the true picture is pro- 
vided by a shaping of the well thereby taking into account the diffuseness 
of the nuclear surface and, in the proton case, allowing the equivalent of 
a variable charge distribution of the nucleus. This has been done by 
W oods and Saxon (1954) for 22 Mev and 18 Mev protons on platinum and 
nickel respectively. A greatly improved agreement was obtained for the 
relation 

V+iw 
IFexpl(r—r)/a] 


VQn= 
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with V40mev, W~ 10 mev, ro(Pt)= 8-24 x 10-8 em, a(Pt)=0-49 x 10-8 em 
and 79(Ni)=5-3 x 10-3 cm and a(Ni)=0:35 X10-%em, This energy region 
allows a more sensitive test of the shaping of the potential well, 
because of the shorter wavelength. It is evident that a multiplication 
of independently variable parameters can be expected to lead to more 
precise agreement between theory and experiment. 
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ABSTRACT 


Two basic equations are derived for deducing the dislocation density in 
powdered materials from the particle size and strain breadth measured 
from the Debye-Scherrer spectrum. In the particle size estimate, 
it is assumed that the material has a block structure similar to that found 
in microbeam studies and that the dislocations lie along the block surfaces. 
The number of dislocations along each face, n, isnot known. In the strain 
broadening estimate the x-ray line broadening from a dislocation array is 
calculated in terms of the broadening due to an isolated dislocation and a 
strain energy factor F', which allows for the effect of dislocation arrange- 
ment. Both methods involve an unknown quantity but by equating the 
two results it is possible in most cases to get both a narrow bracket for 
the dislocation density and considerable information on the dislocation 
arrangement. ‘ 

In annealed metals the values of p range from 2 x 107 cm of dislocation 
line per cm’ for aluminium to 3 x 108 for tungsten and molybdenum, the 
majority of dislocations being ‘ random ’, i.e. the number along each edge 
is approximately 1. There is some correlation between p and the purity 
and annealing treatment of the powder. High purity and high annealing 
temperatures result in smaller dislocation densities. 

In cold-worked metals the dislocation density is more variable. A pre- 
cise analysis is possible for aluminium ; high purity aluminium filed at 
room temperature gives p=4:5 to 7:6 10° with evidence for recovery, 
possibly by polygonization. Commercially pure aluminium filed at both 
room temperature and liquid air give p=5x10° to 2-7x10 and 
p=1-5x 10" to 3:2 101" respectively. Some dislocation array such as 
a ‘ pile-up ’, which causes an increase in the strain energy of dislocations 
exists in both these samples as with all the other metals examined. 
Dislocation densities are also given for tungsten, molybdenum, iron and 
z-brass. In «-brass, the only alloy considered, p is higher than for the 
pure metals and there is strong evidence that the dislocations are effectively 
random and that each dislocation is attached to a stacking fault. 
ee ee eee 

* Communicated by the Authors, 
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$1. INTRODUCTION 


DisLocation densities and information about dislocation arrangements 
have been deduced recently by Gay, Hirsch and Kelly (1958, 1954) from 
X-ray microbeam studies of cold-worked metals. The dislocation 
densities quoted are lower than previously accepted values and it appeared 
desirable to derive values from other x-ray results. In this paper values 
are deduced from x-ray measurements on the Debye-Scherrer spectrum 
of a number of metals and additional information about dislocation 
arrangements is given. 


§ 2. MrerHops 


Two quantities, both of which can be calculated from measurements on 
the Debye-Scherrer spectrum, may be used for calculating the dislocation 
density. The particle size may be derived from measurements of primary 
extinction (Williamson and Smallman 1955), line broadening (Williamson 
and Hall 1953, and Williamson and Smallman 1954) and from measure- 
ments of the peak shift produced by stacking faults (Paterson 1952 and 
Warren and Warekois 1953). The evidence from microbeam experiments 
indicates that the metal is broken up into blocks, the dislocations lying 
in the boundaries between the blocks. If there are n dislocations per 
block face the total length of dislocation line per block is 6 n D/2 (since 
each face is shared by 2 blocks) where D is the dimension of the block. 
The number of blocks per unit volume is 1/D® and thus the dislocation 
density calculated from the particle size is 


eI ad) IL)? eee ee ss s,s, te (2) 


To utilize this equation the value of m must be determined or assumed. 
The value n=1 gives a minimum dislocation density and could apply to 
annealed metals and to very severely deformed metals, in which state the 
dislocations may be very nearly random. 

In deriving eqn. (1) we have defined p as the total length of dislocation 
line in unit volume. With this definition p is independent of the disloca- 
tion distribution. Densities expressed as the number of dislocation lines 
cutting unit area are usually equivalent, but an alternative interpretation 
of this has been employed which is not equivalent to eqn. (1). Gay, 
Hirsch and Kelly (1953), for example, have given values in terms of the 
effective numbers of dislocations cutting unit area with Burgers vector 
lying in the surface of the area. Using this definition one would obtain, 
instead of eqn. (1), that p=n/D?, ie. values } of those quoted in the 
remainder of this paper*. This factor has been allowed for in comparing 
our values with those of Gay, Hirsch and Kelly in that the values quoted 
here are 3 times their published figures. 


* The factor is 4 only for an isotropic distribution of dislocations. For a 
non-isotropic distribution the latter definition is ambiguous and gives a variable 


value for p; by a special choice of the orientation of the area it is sometimes 
possible to make the two definitions agree. 


D2 
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From measurements of line profile it is possible to deduce the breadth 
of the strain distribution ¢ (Williamson and Hall 1953). If S? is the mean 
square strain of this distribution the stored energy of the lattice V is 

V=3 EB S2/2=3 H (Ag*)/2 «ww ee (2) 
where E is Young’s modulus and the factor A depends on the shape of 
the strain distribution. The energy of a dislocation v, can be calculated 
in the absence of interaction with other dislocations and has been shown 
to be, for a screw dislocation 


Le / 
Ue= Ge In (r/r9) toe 8 


where b is the Burgers vector of the dislocation, » is the shear modulus, 
r is the radius of the crystal containing the dislocation and 7 is a suitably 
chosen integration limit, usually about 10-7 cm (see Cottrell 1954). 

In reality dislocations interact so that their strain energies are changed 
on the average by a factor F, ie. 


US peu ke oe 
The dislocation density p calculated from the strain broadening is thus 
V passe 
ae eet y Tl b2 . : . ° ° . “ . (5) 


where k=67HA/u In (r/ro). An approximate value of k can be found 
simply ; H/p is fairly constant for the metals studied here and is approxi- 
mately 2-6, In (7/79) is a slowly varying function of r a reasonable value 
being 4, so that k~12 A. Williamson and Smallman (1954) have in- 
vestigated the effect of line shape and have shown that for cold-worked 
metals the profile of the strain distribution lies between the extremes of 
a Cauchy and Gaussian distribution ; for the Cauchy A~2, and for a 
Gaussian A=47. Thus & lies in the range 2 to 25. 

A more exact calculation of & which does not involve any of the above 
assumptions, is possible by using the results of Wilson (1949) for the 
broadening by a screw dislocation. This calculation is done in the 
Appendix and gives k=16-1 for face-centred cubic materials with Burgers 
vector b along [110], and k=14-4 for body-centred cubic metals with 
b along [111]. These values agree well with the approximate values 
deduced above and support the contention (Williamson and Smallman 
1954) that the actual strain profile is very close to a Cauchy distribution 

Equation (5) can be applied directly to the results only if the interaction 
factor # is known or assumed. The simplest assumption (one that is 
always made in calculating dislocation densities from the calorimetric 
determination of stored energy) is that F=1. One model which 
approximates to #=1 is that in which one dislocation coincides with the 
edge of each block, the dislocations then have their maximum separation 
and interaction is a minimum. Evidently this model coincides with n—1 
in addition to z =1 and it is easy to test the validity of these two assump- 
tions since their correctness would imply that Pp=P,, where Py and p, 
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are the values of p calculated from eqns. (1) and (5) respectively in which 
values n=1 and F=1 have been used. 

This possibility of testing from the experimental results whether a 
given dislocation arrangement is possible can be extended further and 
in addition to giving valuable information about the dislocation arrange- 
ment it considerably reduces the range of uncertainty in the final values 
for p. We shall consider two simple arrangements for which there is 
considerable experimental evidence and which appear to represent the 
two possible extremes for the factor Ff. If other models can be devised 
which produce a greater range of interaction factors these could be used 
equally well in the following discussion. 

Piling up of the dislocations against barriers in the slip plane is com 
monly postulated as a mechanism for producing work hardening and 
this causes a large increase in the strain energy per dislocation. If 
there are n dislocations in each pile up they act to a first approximation 
as one dislocation of Burgers vector nb and hence the strain energy is 
increased from nv, to n?v, (since it is the square of the Burgers vector 
which determines the energy in eqn. (3)). Thus F=n. 

Polygonization on the other hand decreases the strain energy since the 
dislocations ‘ climb ’ from their slip planes to form bend planes and their 
long range stress fields cancel at distances from the bend planes greater 
than h, where h is the separation of dislocations in the bend plane. The 
energy of each dislocation is reduced from the value given in eqn. (3) 
to the value corresponding to r=h. Thus 


F=In (h/ry)/In (7/79). A we ere iO} 


However, h=D/n and substituting eqn. (1) for n, h=3/pD, where p is 
the true dislocation density, also r~}4/3/p and 7y>=10~7 cm.* 
Equation (5) now reduces to an equation for p, namely 


1074/3) / es 
=p n(*) eee LOL ot) meen bee Mies) 
p=p./l=ps O/p 7 

This equation is best solved by trial and error starting with P=Ps 
in the right hand side and, if the results are inconsistent with polygoniza- 
tion (ie. F less than unity) the successive values of p rapidly diverge to 
infinity. If polygonization is shown to be impossible by this treatment a 
lower limit to p can be determined since pile up appears at present to 
represent the extreme degree of interaction. In this case F=n so 
that p=np,=p,/n Thus 


n= (pal pp)? and p=(pspp)””- Sigtema te tt!) 


* This value of , is used in calculating the energy per dislocation to exclude 
the dislocation core for within this region the energy is not proportional to the 
(strain)?. If the value of h closely approaches ry) a value 7)=3 X if on Hy 
be more appropriate since line broadening is sensitive to strain en ae 8 oe 
energy. In the case discussed in later sections h is much greater than 7 anc 


the difficulty does not arise. 
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§ 3. DisLocaTIoN DENSITIES IN ANNEALED METALS 


Annealed metals invariably exhibit extinction. Recent studies by 
Williamson and Smallman (1955 a) show that, for aluminium and a-brass, 
extinction is primary, and will probably be primary in all metals having 
a larger absorption coefficient and similar mosaic size. From measure- 
ments of the fall in intensity due to extinction it is possible to obtain the 
average linear dimension of the mosaic blocks D. 

If a mosaic block is to reflect incoherently with its neighbouring blocks 
it must have a misregister of at least one atomic diameter along each 
edge (see Wilson 1949 a). Exactly such an effect would be produced 
by the block structure assumed in deriving eqn. (1) and the values of 
D deduced from measurements of extinction (or particle size broadening) 
do represent the physical quantity implied by D in eqn. (1). Confirmation 
of this argument, if it were necessary, is given by the very close corre- 
spondence between the values of D found in this study for aluminium and 
the values measured using microbeam methods (Gay, Hirsch and Kelly 
1953). 

Studies of extinction in a number of metals give the values of particle 
D shown in table 1. The value for tungsten has been derived from data 
for the 110 reflection only, the remaining values have been obtained 
using all the observed lines. Values of D for the annealed specimen of 
thoriated tungsten (Warren and McKeenham 1953) have not been included 
in table 1 since D=8x10-*cm giving p,=45x10!; this sample is 
certainly not recrystallized. This result was anticipated because of the 
low annealing temperature used, 1050°c, compared to a temperature of 
1350°c required for purer tungsten in the experiments of Williamson and 
Hall (1953). 

The values for the strain broadening from annealed metals are more 
uncertain than the particle size. Williamson and Hall have shown that 
for aluminium and tungsten the slope of a plot of B* against d* has a 
minimum value (after correcting for the geometrical abberation of 
spectrometer) of about 310-4. The maximum value of this slope is 
probably 4x 10-4. This broadening has been tentatively attributed to 
strains by Williamson and Hall (1953). However, a reliable estimate 
of the strain breadth is difficult since the natural breadth of the K« 
doublet produces an effect indistinguishable from strain broadening 
For copper Kz radiation the half peak width of the wavelength distribu- 
tion is 5x 10-4 A and the apparent strain breadth is 5 x 10-4/1-523 x 10-4, 

Since the maximum experimental error in this determination is 
approximately 1x 10-4, then the probable value for the strain breadth is 
in the region of 0-5x10-4. Taking b as 2:7x10-8cm the value Ps 
(ie. #=1) is about 5x 10-7. This is in very close agreement with the 
values of p, given in table 1 and justifies the appliaction of the block model 
with one dislocation at each block edge in discussing the annealed state. 
Only the values p,, are given in table 1 as these are the more accurate. 

An important correlation also appears between the variations in p and 
the annealing temperature in table 1. This is best seen by comparing 
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the value of p with the ratio 7'y,/7’ where 7'y is the melting point (°K) 
and 7 the annealing temperature. Aluminium, with the lowest value of 
p observed in these studies has the lowest value 7'y/7' (1-21 ) whilst tungsten 
with a temperature ratio (2-28) gives a dislocation density of about 
310-8, The values for molybdenum illustrate this effect most clearly 
since both samples had the same purity and only differed in their annealing 
treatment. This case shows that p is a very steep function of Ty,/T. 
No quantitative correlation using all the results is possible since the purity 
of the material appears to have an appreciable effect. The values for 
iron differ considerably even though both samples had the same annealing 
treatment ; this difference must be ascribed to the difference in carbon 
content since one sample had been decarburized (other traces of impurities 
would be present to the same extent in both specimens). The amount of 
carbon in the armco iron before decarburization impeded the motion of 
dislocations at room temperature to a sufficient extent for a yield point 
to be observed, no such yield point being detected after decarburization. 
However, the magnitude of the effect of the carbon on the value of p 
is at first surprising for at the temperature of annealing the carbon atoms 
would have a very fast rate of diffusion. The effect might be explained 
if it were assumed that, during prolonged annealing after recrystallization 
the dislocations are in constant motion so that the mosaic block pattern 
is constantly changing. 


§ 4. DisLocATION DENSITIES IN COLD-WORKED METALS 


Measurements of strain breadth can be made with comparative ease 
for all cold-worked metals but the values for particle size can be deduced 
with accuracy only for those metals which exhibit extinction even after 
cold working. For the remaining metals the particle size has been 
estimated from the line broadening but this answer is very much less 
accurate. In table 2 two columns each for D are given to indicate the 
reliability. 

The results on aluminium are most useful since accurate particle sizes 
are known. For high purity aluminium the trial and error solution of 
eqn. (7) rapidly converges on a dislocation density p—7-6 10° corre- 
sponding to a strain reduction factor of 1-69, n=162 dislocations per 
surface and to an angle between neighbouring blocks «21°. This tilt 
is in good agreement with microbeam results, as is the particle size D. 
The results are consistent with polygonization having occurred, as 
suggested earlier by Williamson and Hall (1953) from a study of Fourier 
coefficients of the broadened profiles ; but it is interesting that, even in 
this high purity material which recovers readily (so that in some circum- 
stances it recrystallizes at room temperature), the effect of recovery on 
the line broadening is slight, the density is only 70% greater than p,. 
Since the effect of recovery is expected to be less for all the other materials 
studied the values of p, given in table 2 can be taken as upper limits to p. 
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This postulate is confirmed by an analysis of the results for com- 
mercially pure aluminium. Equation (7) has no solution unless h<ro 
and «>20°. Tilts of this magnitude have not been observed in accel 
beam studies and polygonization is very unlikely. Thus Pe= 2-4x 10 
appears to be an upper limit to the dislocation density in this case. The 
lower limit corresponding to severe pile up is p=1-4 x 10° (corresponding 
to n=20). Since commercially pure aluminium recovers some of its 
mechanical properties very readily at room temperature we iis 
anticipate some ‘ unpiling ’ so that p would be in the range 1-2-7 10". 
The results from the aluminium filed at liquid air temperatures are 
consistent with this, although a precise calculation is not possible as D is 
not known. <A crude estimate of D is possible since Gay, Hirsch and Kelly 
(1954) have suggested that D is the slip band spacing, and Brown (1951) has 
shown that for deformation at liquid air temperatures the slip bands have 
one-third the separation they have at room temperature. If Das 
0:7 10-* (Dgp:-=2x 10-4 cm) then there can be no polygonization as 
eqn. (7) has no solution. Equation (8) gives a dislocation density of 
1-5 101° if pile up is severe, whilst the upper limit to the density (p,) 
is 3-2x 10". We expect less ‘ unpiling ’ in this sample and p~2x 101° 
is consistent both with this and the ‘ unpiled ’ value for the sample filed 
at room temperature. 

A similar treatment is possible for the results from tungsten annealed 
at 1150°c which Hall and Williamson (1953) showed not to be recrystallized. 
Solution of eqn. (7) indicates slight recovery, the dislocation density being 
2-6 x 101°, corresponding to a reduction in the strain energy per dislocation 
by 25%, with n=20; and toa tilt « of about 50’ of are. A similar analysis 
could be done using the results from thoriated tungsten published by 
Warren and McKeenhan (1953) since their ‘ annealed’ sample showed 
extinction corresponding to D=8x10-®cm and p,=45x101. It 
seems likely that this sample is polygonized and not truly annealed since 
the annealing temperature was 100°c lower than that used for the recovery 
treatment on the pure tungsten. The strain broadening was not measured 
although the strain breadth must exceed 1:5 10-8. Assuming by analogy 
that nm~10 then p~4-5x 1011. 

It is not possible to carry out such a detailed discussion for the remaining 
results listed in table 2 since exact values of D are not available. It is 
possible to set a lower limit to the particle size by line broadening studies 
and generally 10-*?em >D>10-5em. From the foregoing discussion 
on aluminium we do not expect any polygonization, though even if it 
does occur the effect on the upper limit of the dislocation density will be 
negligible. The range of values of p produced by taking these limits for 
D are given in table 4 for tungsten, molybdenum, iron filed at —183°c 
and iron filed at 20°c respectively. 

Values for p, can also be deduced by the Fourier analysis of line shapes 
which usually gives values of the mean square strain instead of the line 
breadth. The mean square strain for heavily cold-worked metals 
assuming a Gaussian profile is usually one-quarter that obtained assuming 
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a Cauchy profile. The stored energy approach is used to give the two 
extreme values of p listed in table 3. Values for brass are taken from 
Warren and Warekois (1953), who only give values for the Gaussian profile, 
the upper limit is a very approximate estimate from this result. 


Table 3. Calculation of Dislocation Density from the Mean Square 
Strain Determined by Fourier Analysis 


eee ee eee 
| | 


Gaussian Profile Cauchy Profile 
Metal —_———_— ——— eee esate’ 
De Ps S? Ds 
Molybdenum 2x 10-5 So LOEs 8x 105° 1-3 x 1012 
Tron Ie ]0-* I< 1022 G10 1-210” 
a-brass 6x 10-° ie2ec10% 20>¢105° 3x 10" 


The Gaussian and Cauchy profiles bracket the true value of S? and hence 
p; lies between the two values given. 

It should be noted that for «-brass p, is larger than for pure metals. 
However, Warren and Warekois (1953) have shown that there are extensive 
stacking faults in this material and these, in addition to contributing to 
the broadening will act as barriers, maintaining any pile up. For 70/30 
brass they give the separation of stacking faults to be approximately 
120 A and, assuming this to be the particle size, p,—2-5 x 10" suggesting 
that n= 1 and that virtually all dislocations are split into partials separated 
by an extended stacking fault. 


§ 5. Discussion 


For the metals examined in the present study the minimum dislocation 
density in annealed metals varies from 2-7 x 10’ to 3 x 108 and these values 
agree well with the values found using microbeam methods. The 
suggestion (made above) that high temperature anneals produce a lower 
dislocation density is supported by the work of Gay, Hirsch and Kelly 
who quote a value of 1:5 10?cm of dislocation line per cm® for an 
aluminium crystal grown by a method which involves a prolonged anneal 
just below the melting point. The agreement between all these results 
is good; in annealed metals it appears that p can vary from 107 to 
3x 108 em of line per cm*, the value depending on both the annealing 
treatment and the purity. The dislocations in annealed metals are 
arranged in a fairly random fashion, each dislocation forming the edge of 
a mosaic block and producing the tilt necessary to make them reflect 
incoherently. cn 

The probable range of values of dislocation density in cold-worked 
metals is listed in table 4. For heavily cold-worked pure metals the 
values are in good agreement with those of Gay, Hirsch and Kelly and are 
an order of magnitude or more lower than the accepted values. Brass, 
the only alloy considered here appears to have a dislocation density more 
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nearly the conventional value, though this result appears to be associated 
with the presence of faulting. 

In the last column of table 4 the type of dislocation arrangement 1s 
indicated as discussed in the text. The range of values for the dislocation 
density in the tungsten (filed 20°c), molybdenum (filed 20°C) and iron 
(filed 20°c and —183°c) correspond to the uncertainty in particle size 
which has been estimated from line broadening; pin corresponds to 
D=10-4 cm and pyax to D=10-5 em. It should be noted that there is no 
apparent reduction in the dislocation density in tungsten on giving it a 
recovery treatment at 1150°c sufficient to allow polygonization. 


Table 4. Most Probable Range of Dislocation Density in Heavily Cold- 
Worked Metals 


Oe ee 
Dislocation 


Metal Treatment Pmin Pmax Arrangement 


Al 99-99% Filed 20°c | 4.5109 | 7-610°® | Polygonized 


Al (99-7%) Filed 20°c 510° | 2:4101°| Some pile up 
Al (99-7%) Filed —183°o | 1-5x1019| 3-210!!! Some pile up 
W (99-9-++) Filed 20°c | 1-7x1019| 1:7x1011| Some pile up 
W = (99-9+-) | Annealed 1150°c | 2-1 101°) 2-6 101°| Polygonized 
Mo (99-9) Filed 20°c | 15x10] 1-51014| Some pile up 
Fe (99-85) Filed 20°c | 1-05 x 101°] 1-05 x 1011} Some pile up 
Fe (99-85) Filed —183°c | 1-61019| 1-6101!1| Some pile up 
Brass Filed 20°C 2x10" 3x 10!?| Random 
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APPENDIX... 


Calculation of Strain Broadening Produced by a Dislocation in the Powder 
Pattern 


Wilson (1949, 1952) and Frank (1949) have considered the broadening 
produced by a screw dislocation. They find that the broadening is very 
anisotropic, each reciprocal lattice point being broadened in directions 
perpendicular to the Burgers vector of the dislocation but not in parallel 
directions. If the axes of reciprocal space are chosen so that the ¢* 


direction is parallel to the Burgers vector then each spot is broadened to a 
halo whose breadth is given by 


B= ~ (1-41-33 1) 
€9 


where f* is the reciprocal lattice breadth, <) is the effective particle size 
of the crystal cylinder containing the dislocation and h, k and 1 are the 
orders of reflection along the a*, b* and c* axes respectively. The 
component of the broadening due to strains is thus 1-33 J/e, (the breadth 
of the strain distribution is reduced with increase in 7, the radius of the 
cylinder since this increases the volume of material with low strain energy 
hence ¢, appears in the term for the strain broadening). However this 
does not correspond to the broadening observed in practice for two reasons. 
Firstly the experimental broadening is sensitive only to the variation of 
the distance of reciprocal lattice intensity from the origin (i.e. the variation 
in d*). Thus if ¢ is the angle between this line and the c* direction the 
foreshortening factor is cosec ¢ (Wilson 1949). Secondly in a powder 
pattern on the average all possible reflections of each form will 
occur and these will not necessarily have the same value of /. Even in 
hypothetical primitive cubic crystals with dislocation axis along [001] it is 
clear that the 110, 110, 110 and 110 reflections have /=0, so that four 
out of the twelve (110) planes will make no contribution to the strain 
broadening. If ¢ is known, /=d*/b* cos ¢ where b* is the reciprocal of 
the Burgers vector. Thus the effective strain broadening of each reciprocal 
lattice point is : 
B,*=1-33 (d* b*) sin ¢ cos ¢ 
and the strain breadth % (defined as 8,*/d*) of each point is 
1:33 sin ¢ cos d . 
ear cata Rath 16) 

Wilson has calculated the effective particle size for a cylinder of radius, 

r and has shown that 
€g= 167/37 

and hence =0-393 sin 2 ¢. b/r. 7 
In table 5 the various values of ¢, J, and their multiplicities p, 4 and % 
are given for possible reflecting planes from a face centred cubic lattice 
with Burgers vector along the (110) direction. y is the weighted 


mean defined as 4 
p= LXipyp|2p. 
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The mean of the breadth { averaged over lines 3, 4, 11, 12, 16, 19 and 20 is 
€=0-244 bir. 
Since we are considering the case where there is no interaction between 
dislocations D? in eqn. (1) can be equated to zr? and thus 
3 'g 2 
= *(ssize) —=16:127/b?. 


For body centred cubic metals with Burgers vector along the (111) 
direction €=(0-258) b/r and 


oo 8 (ees \ Seep 
Pape) 


Table 5. Calculation of the Strain Broadening Produced by a Dislocation 
in the Powder Pattern from Face-Centred Cubic Metals (Burgers 
Vector along (110)). 


Plane d l p ib b 
111 35° 16 1 4 0:37 b/r 0-185 b/r 
(and 222) 90° 0 4 0 
200 45° 1 4 0-393 b/r 0-262 b/r 
(and 400) 90° 0 2 0 
220 0 2 2 0 0-227 b/r 
60° 1 8 0-34 b/r 
90° 0 2 0 
sii 31° 28 2 8 0-350 b/r 0-268 b/r 
64 46 1 12 0-303 b/r 
90° 0 4 0 
331 ii ae 3 4 0-175 b/r 0-240 b/r 
49° 5] 2 8 0-387 b/r 
Abel l 8 0-241 bir 
90° 0 4 0 
420 18° 26 3 4 0-236 b/r 0-285 b/r 
50° 46 2 8 0-385 b/r . 
718 384 1 12 0-236 b/r 
c= 0-244 b/r 
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IV. Coincidence Measurements on the b-Decay of Boron 12 


By N. W. Tanner 
Cavendish Laboratory, Cambridge* 


[Received August 15, 1955] 


SUMMARY 


A study of B-y and y-y coincidences associated with the f-decay of 2B 
has been made in order to analyse the B-spectrum. It is concluded that, 
relative to the transition to the ground state of 22C, (1:7-0-4)% of the 
f-transitions are to the 4-43 Mev excited state, and <(0-0-+ 0-2)°% to the 
7-66 Mev state. A technique of comparison between !2B and 8Li was 


used which incidentally gave a new upper limit of <(0-2--0-1)% for 
y-radiation following the *Li decay. 


§ 1. INTRODUCTION 

THE f-spectrum from ”B was studied by Hornyak and Lauritsen (1950) 
who showed that 90° of the B-transitions are to the ground state of C. 
For this transition log ft=4:17 which is a normal ‘allowed’ value 
indicating that the 2B ground state is spin 1, even parity ; it is assumed 
that the f-decay comes under the Gamow-—Teller description in accord 
with the isotopic spin selection rule (Christy 1952). Provided J=1+ 
for 12B, then f-transitions to the 4:43 Mev (J =2-+) and 7-66 Mev (believed 
to be J=0-+) excited states of 12C will also be allowed. 

The data of Hornyak and Lauritsen} have been re-analysed in terms of 
the present knowledge of the excited states of *C. An excellent least 
squares fit was obtained with 48° transition to the 7-66 Mev state of 
2C and 3-2% to higher states near 11 Mev ; less than 0-1°% was indicated 
for a branch to the 4:43 mev level. However, plotting N1/?/p(H,—£) 
against # (£ is the total electron energy, Hy the maximum total energy, 
p the momentum, N the counts per unit momentum interval), showed 
that the B-spectrum could, with difficulty, be reconciled with as much as 
a 2% transition to the 4-43 Mev level together with several per cent to a 
higher state, probably that at 9-61 Mev. 

Coincidence measurements on the 8 and y radiations from ’B have 
previously been made by Vendryes (1951). Unfortunately the work was 
done with Geiger counters which did not allow a distinction of B~y 
coincidences from £-bremsstrahlung coincidences with any certainty. 
Vendryes’ conclusion was that (4+1)% of all the B-decays are to a 
y-decaying state of ?C, probably that at 4-43 Mev. 


*Communicated by the Author. 
+I am grateful to Professor T. Lauritsen for the opportunity of studying 
the original data on the f-spectrum, 
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§ 2, Propuction or ?B 

The convenient manner of producing !2B is by the reaction 4B(d, p)B : 
a thick natural boron target gave adequate yield with the deuteron energy 
and beam available from the Cavendish 1 Mv accelerator. 

As the half-life of 12B is short (25 m sec), and the background from a 
deuteron bombardment very intense, it was considered essential to use 
a pulsed beam technique, ie. alternate bombardment of target and 
counting of radiation. An elementary calculation showed that maximum 
efficiency of !2B counting would be achieved for a pulsing cycle extending 
over about two half-lives (50 m sec) or less. 

The arrangement used consisted of a cylinder mounted within the 
vacuum system in the path of the deuteron beam. A synchronous motor 
rotated the cylinder, via a Wilson seal shaft, at 300 r.p.m. Four slots 
in the cylinder allowed the beam to bombard the target during alternate 
25 msec periods. A lamp and photocell arrangement, using the same 
rotating cylinder as shutter, provided a signal to the master timing circuit 
which allowed counts to be recorded only while the beam was off the 
target. 

The master timing circuit consisted essentially of electronic gates which 
intercepted pulses in appropriate parts of the counting circuits (see fig. 2). 
An additional function (not shown in fig. 2) performed by the master 
timing circuit was to switch off the scintillation counters while the 
deuteron beam was bombarding the target. The immense count rate 
during bombardment was liable to cause instability in the photomultiplier 
gain. Switching of the photomultiplier was obtained by applying a 
120 volt, 20 msec long rectangular pulse to the first dynode; during 
the pulse the voltage of the first dynode was pulled below that of the 
second dynode so that the photomultiplier gain was zero. 

In the following sections, no further mention will be made of the pulsed 
beam technique, as this system was used for all measurements. 


§ 3. COMPARISON OF 2B AnD SLi 


Any attempt to detect y-radiation following the B-decay of !2B, or any 
other high energy 6 emitter, has special difficulties. Absorbers can be 
used to prevent the y-counters detecting the B-particles directly, but only 
at the cost of a considerable intensity of bremsstrahlung (even for materials 
such as wax or graphite). Anti-coincidence techniques are possible, but 
in practice it was found difficult to achieve the necessary 100° counting 
efficiency for the £-detector. 

A unique solution to this difficulty was found in the remarkable 
similarity of the B-spectrum of SLi to that of !B, together with the 
knowledge that there is little or no nuclear y-radiation associated with 
the *Li decay. Bunbury (1953) has put an upper limit of 0-8°% on the 
existence of B-y coincidences from 8Li and this was confirmed in the 
experiments to be described, 
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Figure 1 compares the two £-spectra normalized to the same total 
count. The ”B spectrum has rather more high energy £-particles than 
the *Li spectrum. The agreement, however, is sufficiently good to justify 
using *Li for measuring the bremsstrahlung background while looking for 
nuclear y-radiation from 12B. An approximate calculation, described in 
the Appendix, showed that the !2B f-particles would give a bremsstrahlung 
intensity per 6-count in excess of that from ‘Li by about 21°. 


Fig. 1 
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Comparison of the f-spectra of 2B and 8Li normalized to the same total count 
for momenta greater than the discriminator cut-off. (After Hornyak 
and Lauritsen 1950.) 


Thick Li,SO, targets were used for producing the 0:9 sec *Li activity 
via the “Li(d, p)§Li reaction. A target mounting was arranged, with 
lithium on one side and boron on the other, which could be reversed without 
opening the vacuum system. For a given deuteron energy and SE 
the equilibrium f-activities of *Li and ”B were very much. the same. 
The targets could be distinguished by the lifetimes of the activities. 


§ 4, B-y COINCIDENCES 
The arrangement of electronics and counters is illustrated in fig. 2. 
Detectors used were sodium iodide scintillation counters and a small, thin- 
walled proportional counter filled with argon to 60 cm pressure. A Hees: 
of especial interest was the use of a proportional Pe ee 
counter telescope for detecting the energetic B-particles. lue scintil een 
counter allowed the use of a fast coincidence circuit (0-02 usec) with the 
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invaluable reduction of random coincidences. At the same time, the 
proportional counter distinguished genuine f-particles from y-radiation 
(bremsstrahlung) or neutron-induced 181 activity in the Nal crystal 
itself. 

Pulses from the y-scintillation counter were gated to a Hutchinson- 
Scarrott pulse height analyser, the gating condition being that simultaneous 
counts (within the resolving times indicated) are recorded in the propor- 
tional counter and f-scintillation counter. Use of the pulse height analyser 
is vital to the experiment, as it is only in this way that the energy of the 
y-radiation can be identified. Energy calibration was obtained using 
ThC” and Po—Be sources. In addition the Po—Be was used for absolute 
efficiency calibrations as its activity was known to a 10%, accuracy. 
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Block schematic diagram of electronic and counter arrangement for By 
coincidence measurements. 


Counting was carried out one hour at a time, alternate hours being used 
for a total coincidence count and a random coincidence count ; the latter 
was arranged by the introduction of a 0-05 psec delay one side of the fast 
coincidence circuit. Coincidences between the proportional counter and 
B-scintillatiom counter were recorded on the scaler and served to monitor 
the B-activity. The spectrum on the pulse height analyser was recorded 
photographically. For those runs in which comparison with 8Li was made. 
the ®Li, 12B counts were alternated. 

Initially measurements were made with pulse height discrimination 
such that only those f-particles with energy >2 ev were recorded. 
Even without SLi for comparison, a 4:4 Mey y-line could be identified. 
Figure 3 shows coincidence y-spectra from ”B and 8Li, normalized to 
the same f-count, together with the calibration y-line from Po-Be. The 
apparent 4-4 Mev y-radiation could arise by a f§-transition either to the 
4:43 Mev level of 12C or to the level at 7-66 Mev with subsequent y~y 
cascade decay. In the latter case a 3-2 Mev y-radiation must also occur. 
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From the data shown in fig. 3 an upper limit could be placed on the 
occurrence of a 3-2 Mey radiation and hence a 6-transition to the 7-66 Mev 
level of ’C. Relative to the ground state f-transition this limit is =0-39%, 
Assuming zero transition to the 7-66 mev level, a (1-:7-+-0-4)°% transition 
to the 4:43 Mev level is indicated. These figures represent the mean of 
several measurements interpreted in the light of the 8Li/!2B comparison 
spectra of fig. 3. Correction was made for the effect of the finite dis- 
criminator bias on the relative counting efficiency for ground state and 
excited state B-spectra. 

Figure 3 also serves to justify the use of *Li for comparison, in that the 
spectrum indicates that <(0-2-++0-1)°/ of the SLi B-transitions are followed 
by a ~4-9 Mey y-radiation. 
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Spectrum of y-radiation in coincidence with f-particles of energy greater than 
2mev. (Negative counts are a result of subtracting random coincidences.) 


Further By coincidence spectra were taken with the discriminators set 
to exclude f-particles with energy <5-7 Mev, i.e. the end point energy of 
a 12B B-transition to the 7-66 Mev level of "C. Statistics of these counts 
were so bad that they added little but qualitative evidence to the 
impression that a B-transition does take place to the 4:43 Mev state of AC: 

A simple re-arrangement of the electronics served to display she 
B-spectrum, from the f-scintillation counter, on the pulse height analyser. 
However, even with the coincidence condition of a y-radiation of energy 
>3-3 Mev, no useful information was obtained. The broad nature of he 
B-spectra makes the identification of end point energies very inexact 
unless good statistical accuracy can be achieved. 


§ 5. y-y COINCIDENCES 
The technique for this experiment involved the use of two scintillation 
counters (both with graphite f-absorbers), fast coincidence circuit and 
pulse height analyser, similar to that described in §4. In this case, 
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however, the output from one counter was fed to the analyser under the 
conditions that a coincidence occurred and that the pulse height from the 
other counter corresponded to a y-ray of energy between 3-3 and 4:5 Mev. 
The proportional counter was used to monitor the B-activity. 

The problems of bremsstrahlung background and low coincidence rates 
are very much greater in the case of yy as compared to fy measurements. 
Even with the counters set 180° apart (an obvious minimum for 
bremsstrahlung effects and the expected maximum for y-y coincidences 
from the 7-66 Mev state of 12C (Seed 1955)) the bremsstrahlung coincidences 
were sufficient to mask completely any nuclear y-radiations. This was 
clearly demonstrated by a measurement of the angular correlation which 
fitted the form cosec? 6/2, i.e. coincidence rate proportional to the inverse 
square of distance between the two counters: perhaps explicable in terms of 
Compton scattering of a bremsstrahlung photon or double bremsstrahlung 
interaction by a single electron ; nuclear y—y coincidences are quite clearly 
excluded. 


Fig. 4 
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Spectrum of y-radiation in coincidence with another y-radiation of energy 
between 3:3 and 4:5 Mev. (Negative counts are a result of subtracting 
random coincidences.) 


Rather more interesting results were obtained by the 8Li/!2B comparison. 
Figure 4 shows the coincidence y-spectra normalized to the same monitor 
f-count. The main characteristic of the two curves is their similarity 
although the B count is significantly greater. Taking the difference 
between the two curves gives an upper limit value for a 2B(Byy)!2C process 
relative to 1B(8)"C as 0:-47%. If however we allow 21° more 
bremsstrahlung photons per B-count from 12B as compared to SLi (see 
Appendix), then the upper limit for a nuclear y-y cascade is reduced to 
(0-04 0-20), where the error is that from counting statistics alone. 
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§ 6. Discussion 


The evidence presented indicates that the B-decay of 2B goes primarily 
to the ground state of #2C but with (1-7-+-0-4)°% branch to the 4:43 Mev 
level. Transitions to the 7-66 mev level are excluded to the order of an 
upper limit of (0-04+-0-20)%, provided that the state does decay only by 
y-y cascade. Uebergang (1954), Beghian ef al. (1953) and Seed (1955) all 
support this contention (the factor of three discrepancy between Uebergang 
and Beghian et al. is merely instrumental in that it arises from the use of 
a thick target by the latter workers). Rasmussen ef al. (1955) on the 
other hand conclude from a study of inelastic scattering of «-particles 
by #C, that the 7-66 mev level of C decays to «-particles with a probability 
of at least 80% ; however, it is not clear that their observations might not 
have been confused by a strong angular distribution and an unfortunate 
choice of angle. 

The allowed f-transition to the 4:43 Mev state of !2C (allowed since !2B 
ground state is (1, +) and "C first excited state is (2, +-)) has log ff=5-1. 
Assuming the 7-66 Mey level of 12C decays solely by y—y cascade, then the 
B-transition to this state has log ft>5-6. 

The interpretation of the f-spectrum (Hornyak and Lauritsen 1950) 
has already been discussed in § 1. From the coincidence results, it seems 
- that analysis in terms of a 5% transition to the 7-66 Mev level of ?C, is 
quite definitely incorrect. 
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APPENDIX 
Difference of Bremsstrahlung from *Li and ”B 


Heitler (1947, p. 168) gives the extreme relativistic approximation to 
the cross section ¢(k, #) for the emission of a bremsstrahlung photon 
of energy between k and k+-dk by an electron of energy H as 


g(t, By~[ A +5 | [toe 3]: 


We are interested in a particular value of k, say 3 Mev, but all values of # 
from initial energy H=E, down to # =k as the electron loses energy by 
inelastic collisions in the absorber (ionization loss is far more important 
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than bremsstrahlung loss for ~10 mey electrons in graphite). Hence the 
required cross section is 


r Eo 
b(k, EB) dE. 


@(k, Ey)= I 


Unfortunately the integral cannot be handled analytically, but a graphical 
method gave the result 
O~(H,—k)**. 


It is now necessary to integrate the product of @®(k, Hj) and the 
function N(E,) dE, which describes the B-spectrum of ‘Li or *B shown 
in fig. 1. The integrals obtained, again by the graphical method, are 
proportional to the number of bremsstrahlung photons of energy k& per 
B-particle count. The ratio of the 12B integral to that for SLi was 1:21 
for k=3 ev. 
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SUMMARY 

A description of a practical interferometer system based on the principle 
of post-detector correlation is followed by a short account of an instrument 
incorporating this principle. This instrument has been used at Jodrell 
Bank for measurements of the angular diameters and structures of the 
radiation from the quiet sun and the sources Cygnus and Cassiopeia. 
The method of computing the results is indicated, and a short account 
is given of a synthetic source which is used to calibrate the equipment. 
Reference is made to a method of obtaining the gradient of the correlation 
function by the use of three adjacent operating frequencies. The results 
of the measurements on Cassiopeia and Cygnus are given in Part II of 
_ this paper. 


$1. InrRopDUCTION 


Ear_y attempts to measure the angular diameters of the localized radio 
sources (radio stars) by using interferometer techniques were unsuccessful, 
and only an upper limit of six minutes of arc could be placed on their size 
(Ryle and Smith 1948). In anticipation that the sources might subtend 
angular diameters comparable to those of the visible stars, a new inter- 
ferometer technique was proposed by Hanbury Brown for measuring 
the angular diameter of the two most intense radio stars, Cygnus and 
Cassiopeia. The present paper describes the development of the inter- 
ferometer and of a small prototype of the instrument which was used to 
test the principle in a model experiment to determine the diameter of the 
solar disc at a frequency of 125 Mc/s. An improved full scale instrument was 
completed in 1952 and applied to the measurement of the angular diameters, 
and later the fine structures, of the sources in Cygnus and Cassiopeia. The 
results of these measurements are reported in part IT of this paper. 

The principle and the theory of the interferometer have been discussed 
in detail by Hanbury Brown and Twiss (1954). The system does not 
require coherent local oscillators in the two receivers but relies upon the 
coherence of the noise spectra formed after the signals in each receiver 
have been fed into separate square law detectors ; for this reason it has 
been termed a post-detector interferometer. The instrument may be 
used on very extended baselines ; it is capable of a much higher resolving 
power than previous interferometers and is not seriously affected by 


ionospheric scintillation. 
* Communicated by Professor A. C. B. Lovell, F.R.S. 
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§ 2. THe Post-DeTecTor INTERFEROMETER 


The basic arrangement of the post-detector interferometer is illustrated 
in the schematic diagram of fig. 1. The aerial systems A, and Ag, spaced 
at a distance d apart, are connected to two independent receivers having 
identical characteristics. The receivers are double superheterodynes, 
both tuned to 125 Me/s and have a radio frequency bandwidth of 200 ke/s. 
The bandwidth is accurately determined in the final I.F. strip of each 
receiver immediately prior to a square-law detector. 
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Block diagram of the post-detector interferometer. 


The application of a ‘ white’ noise spectrum to one of the non-linear 
detector elements produces a system of beat frequencies in the output. 
These beat frequencies are then passed through a low frequency filter 
of pass band 1 ke/s to 2-5 ke/s. The outputs from both receivers are 
treated in a similar manner and the two low frequency outputs are then 
fed into a common mixer or correlator unit. The function of the correlator 


is essentially to record the r.m.s. levels, (Z,?)? and (#,?)1/2, and the 
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scalar product, (H,E,), of the two sets of signals present in the low 
frequency outputs of the two receivers. From these three parameters 
the cross correlation coefficient (p2) of the two low frequency spectra may 


be computed. It can be shown (Hanbury Brown and Twiss 1954) that 
the function 


te Ff? sin (w/c) + F? cos (wd /c) 
FF? cos (0) 
where w, is the angular frequency to which the aerials are tuned, d is 
the baseline of the interferometer, and c is the velocity of light. The 
functions ’;, and F,,, refer to the sine and cosine terms of the Fourier 
transform of the angular distribution of intensity across the source. 
Hence by varying the spacing of the aerials, the angular diameter of 
the source may be computed, and subject to assumptions of symmetry, 
the intensity distribution may be determined. 

In operation the two aerials A, and A, are each fixed in azimuth and 
in the same plane of polarization and are adjusted in elevation to the 
declination of the source. As the source transits across the aerial beam, 
the increments in the root mean square values of the low frequency outputs 
together with their scalar product are recorded on a three pen recorder. 
The value of the cross-correlation coefficient is then obtained from an 
analysis of the three records. 

When a large aerial spacing is used the low frequency output of the 
distant receiver is transmitted over a radio link by modulating a high 
frequency carrier. The time of arrival of the two low frequency outputs 
at the correlator is equalized by inserting in the output of the local 
receiver a delay equal to the time of transmission along the baseline. As 
all the parameters in the receivers are identical, the recorded correlator 
output of the product will be a maximum when the signals entering the 
receiving aerials are identical, but will diminish with the intensity of 
either and with the phase shift between their spectra after detection. 
It is essential that all the components in the equipment should preserve 
the relative phase of the two low frequency signals, but the problem is 
considerably simpler in practice than the preservation of the relative 
phase of the radio frequency signals in conventional interferometers. 

The background level of noise in the two receivers is entirely uncorrelated 
and hence does not contribute to the mean value of the correlation. The 
effect is, however, to produce a random background of noise in the 


channel recording (H,H,), and it is largely this factor which finally limits 
the sensitivity and resolution of the equipment. 


§ 3. DESCRIPTION OF THE APPARATUS 
(a) Aerial System 
A high level of resolution may be obtained with the post-detector 
correlation interferometer provided that the signal to noise ratio in each 
channel is of the order of unity, but the resolution falls off severely at lower 
signal levels. In order to achieve a signal to noise ratio of the order of 
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unity from the radio sources Cassiopeia and Cygnus on an operating 
frequency of 125 Me/s each aerial was required to have a gain of 600 over 
a half-wave dipole. 

The necessity of moving one of the aerial systems to various sites, In 
order to determine the variation of correlation with baseline separation, 
dictated the adoption of a prefabricated antenna which could be readily 
dismantled and re-erected. The most satisfactory solution was found to 
lie in a collinear broadside array of full wave dipoles, comprising 320 
separate elements fed in 16 bays by equal lengths of coaxial cable. The 
overall dimensions of each aerial were 128 ft. x 42 ft. and the corresponding 
beam widths between half-power points were 24° in Right Ascension and 
11° in declination. The aerials could be tilted about central tripods so 
that both Cassiopeia and Cygnus were received in the centres of the 
beams. When it was necessary to dismantle one of the aerials and re-erect 
it on a new site the whole structure could be collapsed into 32 separate 
flat units which were readily transportable. 


(b) Receivers 

The main requirements for the receiving systems were that the band- 
widths and local oscillator frequencies should be accurately matched. 
The law of detection had to be the same in both channels and the law 
itself had to be accurately known or calibrated. In order to maintain 
equal bandwidths in the two receivers a double superheterodyne 
system was employed in which the gain was accomplished in a broad 
band intermediate frequency amplifier at 10 Me/s, whilst the bandwidth 
was dictated by a single damped tuned circuit at 2-1 Me/s. These tuned 
circuits were accurately matched in the two receivers and were arranged 
to give identical bandwidths over the whole of their response down to 
—40dB. The overall bandwidth at —3 dB was 200 ke/s. 

The square law detectors each consisted of two pentode valves, type 
6B8, operated with the screens regulated at 105 volts and the control 
grids biased at —13-5 volts. The signal was applied in antiphase to 
the two control grids whilst the output was taken across a common load 
of 10K in the anode circuit. This arrangement rejected the first, 
third and fifth order harmonics whilst the low level of the fourth order 
yielded an accuracy of approximately +0-5% over an input excursion 
of +5 volts. 

The final stages of the receivers were filter units, defining a low 
frequency bandwidth of 100-6000 c.p.s. from the output of the square law 
detectors, and associated amplifiers and cathode followers to feed in one 
case into the relay transmitter and in the other case directly into the 
correlator. 

(c) Relay System 

The output of the receiver at the distant station consisted of signal and 
noise frequencies in a band from 100-6000 c.p.s. These were relayed to 
the home station by a radio link operating on a frequency of 83 Me/s. 
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The relay transmitter was crystal controlled and could be operated at a 
maximum input of 300 watts. The slight phase dispersion occurring 
in the modulator was compensated by a network having the same 
characteristics in the opposite limb of the equipment. The radio 
frequency output of the transmitter was fed into a vertically polarized 
corner aerial which was beamed towards the home station. 

The relay receiver was a superheterodyne, crystal controlled to the 
transmitter frequency and equipped with a powerful automatic gain 
control. 


(d) Compensation Network and Final Filters 


It was essential that the signals fed into the final correlator unit should 
be accurately contained in identical bandwidths and that the phase and 
amplitude dispersion over these bandwidths should be a minimum. The 
relay system introduced a slight dispersion in one channel and this was 
compensated by a network of identical characteristics in the opposite 
channel, where provision was also incorporated for introducing an 
adjustable time delay to compensate for the delay in the transmission of 
the signal between the two stations. The delay lines were tapped ladder 
networks working into matched loads and were capable of a maximum 
delay of 300 see with a level response up to 10000 c.p.s.; they were 
originally designed for use with anticipated baselines of the order of a 
hundred kilometres. The delay incurred with a baseline of 12 km., the 
longest baseline involved in the present experiments, was 40 sec 
corresponding to a phase shift of 14° at 2000 c.p.s. 

At this stage the final video bandwidth was defined between 1000 
and 2000 c.p.s. by accurately matched resistance—capacity ladder filters. 
An improvement in the general signal to noise ratio could have been 
achieved by increasing this bandwidth, but only at the expense of further 
complications incurred in maintaining a low phase dispersion. 


(e) The Correlator 

The correlator accepted the signals from the final filter units in the 
band from 1000 to 2500 ¢.p.s. and converted them into their root mean 
square readings for display on two of the pens of a three pen recorder. 
To the third pen it fed a signal corresponding to the mean value of the 
product between the two inputs. The r.m.s. recordings were obtained 
by rectifying the signals and integrating with a suitable time constant. 
The rectified signal level was arranged to be sufficient to give full scale 
deflection of the recording instruments without the need for d.c. 
amplification. 

The operation of the multiplier unit may be understood by reference 
to the diagram in fig. 2. The outputs of the final filter units were con- 
nected to the two points A and B, whence one signal was caused to suffer 
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a reversal of phase every 1/30th second before they were combined in an 
additional mixer. The combined signal at the point C therefore 
alternately consisted of the sum and difference of the two spectra at 
A and B, 

(£3) + (Ey), JL) s+ (Lo)n] : . - 2 (1) 
where (H,), and (H,), refer to the low frequency signals and (#,), and 
(E,), to the corresponding noise terms in the two original channels. 
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The cross multiplier circuit of the correlator unit. 


Upon squaring the signal expressed in eqn. (1) in the following push—push 
square law detector, and after filtering out the high frequency terms, 
there remained a constant term proportional to the product 

(HY) < (Bs), . a 
the polarity of which reversed every 1/30th second in synchronism with 
the switching waveform. The remaining terms all produced fluctuating 
signals or were constant without variation at the switch frequency and 
were largely filtered out in a selective amplifier before being finally 
presented to a phase sensitive detector. The effect of the latter was 
to take the difference of the two terms expressed in (2) and, after integra- 
tion in a simple filter, to feed to the recording meter a signal directly 
proportional to (#,), x (#.) 


st 


§4. ASSESSMENT OF CORRELATION, SYNTHETIC Norsk SouRCE 


The appearance of the chart after the transit of a source is depicted 
in fig. 4. The lower record shows the cross modulation product between 
the components forming the upper two records. As the source is resolved 
the amplitude of the lower record is reduced whilst those of the upper 
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records remain unchanged. The correlation coefficient (p?) may be 


computed by measuring the amplitude of the increments (#7)? and 
(Z,?),’’? in the two receiver channels and comparing their product 
with the observed increment from the third output of the correlator. 

Assessment of correlation in this manner requires that the constants 
of the equipment be accurately known and invariable over long periods. 
A more reliable method of assessment, which also avoids any possible 
error in the whole of the correlator system and final filter units, is to 
introduce, immediately upon the completion of the transit and under 
identical conditions, a synthetic noise signal whose correlation is known 
to be unity. 
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The synthetic source. Adjustments of signal and background noise are made 
on the controls R,, R., R; and R, to the left of the diagram. 
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The correlated component of the synthetic noise signal was produced 
at a frequency of 200 kc/s in a neon stabiliser tube, whence it was fed 
separately via two potentiometers, R, and R, in fig. 3, to two indepen- 
dent but identical radio frequency filter units, having a passband of 
150 ke/s at 300 ke/s. Uncorrelated noise signals from the other two 
neon generators were added independently to each channel through the two 
gain controls R; and R, prior to passing through the filter units. The 
envelope of the noise in each channel was detected by two square law 
detectors similar to those in the main receivers, and the output frequencies 
were limited to a similar band by two audio frequency filters of identical 
characteristics. 

The two outputs of the unit were applied to the correlator immediately 
after the transit of the radio star and the four controls R,, Ry, R3, and 
R, were adjusted to give identical conditions in each channel to the 
actual signal and background noise occurring at transit. An example 
of the application of this method of calibration is shown to the right of 
the transit of the source in fig. 4, which depicts a source of 0-65 correlation 
compared with a synthetic noise signal of unity correlation. The ampli- 
tude of the rise in the two intensity recorders is the same for the true source 
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and the synthetic source, but the cross modulation product in the latter 
case is of almost double the amplitude. 


Fig. 4 
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A typical record of the transit of the source in Cassiopeia. The synthetic 
calibration to the right of the transit shows that the Cassiopeia source 
had a correlation coefficient p?—0-65 at this aerial separation. 


The operation of the synthetic source itself was checked by comparing 
the correlated output with the correlation of noise in a single audio 
channel fed into both inputs to the correlator. This enabled the unity 
point to be compared; the zero point was checked by applying large 
uncorrelated inputs in the absence of a correlated signal. 


§5. FrREQuENCcY CHANGING 


The resolution of the interferometer depended upon the number of 
wavelengths separating the two aerials. As the aerials were moved 
apart the resolution increased and the correlation coefficient obtained from 
the records varied according to the square of the Fourier transform of the 
intensity distribution across the source. Instead of varying the separa- 
tion of the two aerials a similar change could be achieved by maintaining 
a constant aerial separation but varying the frequency of the radio noise 
selected by the interferometer from the radiation of the source. 

Provision was incorporated in the receivers for simple frequency changing 
so that they could be operated on 116-5, 125, or 132-5 > Me/s, as desired. 
The three values of the effective baseline (measured in wavelengths) 
so obtained differed by about 12° and afforded a useful means of 
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measuring the gradient of the correlation function without the need to 
move the remote aerial system and its associated equipment. 

In utilizing this principle it was assumed that the intensity distribution 
of the radiation across the source was not sharply dependent upon 
frequency. This assumption was justified for the source in Cassiopeia 
and Cygnus by noting the effect on a gradient of the transform which had 
been determined previously by variation of the aerial spacing. Further 
confirmation is also available by comparing the results of Smith (1952) 
and Mills (1952) on frequencies of 210 and 101 Mc/s with those obtained 
on a frequency of 125 Mc/s with the present interferometer. 


§6. THe EQUIPMENT IN OPERATION 


A prototype of the interferometer was completed in 1951 and the 
principle of its operation was fully tested with aerial apertures of 
100 sq. metres by measurements of the distribution of radiation across 
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the solar disc on a frequency of 125 Mc/s. The measurements were made 
at the time of considerable sunspot activity and this contributed a large 
random error in the readings, together with a systematic increase in the 
level of correlation at large aerial spacings. A plot of the observed corre- 
lation at a number of aerial spacings is shown in fig. 5. The results 
obtained were consistent with the distribution determined by more 
conventional techniques, and after minor modification and increase of the 
aerial aperture the equipment was applied to the measurement of the 
diameters of the cosmic radio sources Cygnus and _ Cassiopeia. 
An extensive series of observations was conducted on these two sources 
and a detailed account of these observations is given in Part II of this 
paper. It was found that the equipment operated satisfactorily under 
all conditions of the experiment and its accuracy was not seriously affected 
by the presence of severe scintillations. Since both Cassiopeia and 
Cygnus were resolved at distances much shorter than the extreme base- 
lines which could have been used with this interferometer, the full 
potentialities of the technique have not yet been realized. 
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SUMMARY 
An account is given of the application of the post-detector interferometer 
to measurements of the radio stars Cygnus A and Cassiopeia A. An outline 
of the experimental procedure and results is followed by a brief comparison 


of the photographic and radio frequency distribution of intensity of each 
source. 


§ 1. LyTRopUCTION 

THE post-detector interferometer described in the first part of this paper 
was applied to the measurement of the approximate angular diameters of 
the intense radio sources in Cygnus and Cassiopeia, and upon the successful 
completion of this programme it was used to determine the angular 
distribution of intensity across these sources. 

A brief review of the relevant existing information on Cassiopeia and 
-ygnus is given in table 1. 


§ 2. MEASUREMENTS WITH THE Post-DETECTOR INTERFEROMETER 
(a) Initial Testing and Calibration of the Instrument 


The new technique was initially tested by performing a model experiment 
to determine the angular diameter of the solar disc (see Part I). The 
improved full-scale version of the instrument was then set up on a baseline 
of 300 metres north-south in order to check the performance of the 
equipment when recording Cassiopeia and Cygnus. 


(b) The Operational Method 

All the measurements described in the succeeding sections were con- 
ducted in the following manner. The aerials were adjusted to the 
declination of the Cygnus source, which was allowed to transit. The 
receivers were then disconnected and the correlator was calibrated by 
means of a synthetic source. Two uncorrelated signals were applied and 
adjusted to give the same background noise as those observed. After a 
Ge ee ieee Sa Ee 
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Table 1 
Cygnus Cassiopeia 
Flux* 
watts M-? (c/s) stipe el Mei 2:2 10-*? 
at 81 Me/s 
Radio 
source Positiont 19b 5]™ 458+ 18 23h 21m 128+ 18 
R.A. 
Dec. 40° 35’+1’ 58° 32-1'+0-7' 
Filamentary 
Identified Galactic emission 
object collision nebulosity 
Photographic Photographic 15-05 — 
identification? magnitude mpg 
Angular diameter 18” x 30” 6’ 
Distance 3:3 107 parsecs | 200-2000 parseecs 


(old scale) 
* Ryle e¢ al. (1950). + Smith (1951). { Baade and Minkowski (1954). 


short time a correlated signal was also applied and adjusted so that when 
observed on the two r.m.s. intensity recorders the level on both channels 
coincided with that observed from the source. The resulting level on the 
cross multiplier was then measured and the whole calibration procedure 
was repeated. 

The source in Cygnus presented a complication in the uncertainty of 
the relative levels of the intensity of the radiation from the source and 
that from the Galaxy. The source is situated on a steeply inclined gradient 
of the galactic radiation, the uncertainty of the true contour of which 
masks the true level of the base of the Cygnus transit. The correct level 
was eventually established by accurately calibrating the contour of the 
narrow aerial beam in R.A., and also by a detailed analysis of the cross 
multiplier records, from which the contributions from the localized source 
could be assessed at all parts of the transit. 

A similar procedure to that used on the Cygnus source was carried out 
for the source in Cassiopeia. Between ten and twenty runs were made 
on both of the sources at each site and every run was calibrated by means 
of the synthetic source. Throughout the series of measurements the 
receiver and relay transmitter were frequently removed to the ‘ home ’ 
station at Jodrell Bank and tested. Apart from normal testing with 
signal generators, a method of testing was developed which immediately 
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showed up any loss of correlation in the whole equipment. This method 
was to split the home aerial into two halves and to apply the two receivers 
to these. The transit of the Cassiopeia source was then recorded and the 
equipment was calibrated with the synthetic source to check that it showed 
unity correlation. The synthetic source itself was checked by reducing 
the uncorrelated component to zero and comparing the resultant correlation 
with that obtained by connecting together both the inputs of the correlator 
and feeding it with the same level of noise from asingle channel. A further 
check was made by removing the correlated signal and increasing the two 
uncorrelated background noise levels, so that in the absence of a correlated 
signal the two r.m.s. intensity recorders registered a high reading, whilst 
the cross multiplier was not deflected. 


(c) The Preliminary Diameter Measurements 

The diameter measurements were made on a north-south baseline of 
3-9 km with the remote station in a direction 1° east of south. The results 
showed that the radio sources were both partially resolved. The corre- 
lation of the Cygnus source had fallen to about 0-8 whilst that of the 
source in Cassiopeia had fallen almost to zero. A number of readings 
were taken to ensure that the results were consistent and the angular 
widths of both the sources, expressed as the angular widths of equivalent 
strip sources of constant surface intensity were_then calculated. These 
results, which have already been published (Hanbury Brown et al. 1952), 
showed that the source in Cygnus had an angular diameter of the order 
of thirty-four seconds of arc, whilst that in Cassiopeia had an angular 
diameter greater than two minutes of arc. 


(d) The Test for Radial Symmetry 

The initial measurements of approximate angular diameter were 
followed by an investigation into the general shape of the sources. The 
next series of measurements, which were conducted in the autumn of 1952, 
were therefore made at a constant aerial spacing along three different 
azimuths, one of which was chosen to coincide with the north-south 
azimuth used in the initial diameter measurements. The separation of 
the aerials was decreased from 3-7 km to 2-16 km in three azimuths 
inclined to each other at angles of approximately 120°. The reduced 
baselines were used on the assumption that it would then be possible 
to obtain correlation coefficients for both sources which were on the steep 
gradient of the first maximum of the function and which would therefore 
permit a sensitive measurement of changes in the primary distribution 
with azimuth. 

The values of the correlation coefficients and the equivalent angular 
strip diameters determined by the measurements along the three axes 
are given in table 2. 
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Table 2 
Aerial spacing Cygnus Cassiopeia 
Bearing : : Angular Cieeint Angular 
Distance | degrees Cor ce width of oe eer width of 
in km from N. SEU equivalent | °°° eae equivalent 
(clockwise) pP strip - strip 
2-16 113 0-30+0-03 | 2’ 10"+4”" | 0-:08+0-02 | 2’ 55" +10" 
2-16 235-5 0-79+0-08 | 1’10"+7" <0-01 +3’ 30” 
2-16 176 0-97 + 0-04 <o: 0:05+0-03 | 3’ 20"+15" 
3-99 177 0:79+0-07 | 0’ 34”+8" <0-07 unreliable 
0-30 349-5 0-99 + 0-10 cad 0:96+0-09 | 4'30"+1’ 


* The value of p? is defined in Part I of the present paper. 


From these readings it was at once apparent that the source in Cassiopeia 
was a roughly symmetrical object whilst that in Cygnus was markedly 
asymmetrical. The Cygnus source was considerably more extended in 
a direction approximately east-west but it was found that no simple 
elliptical distribution would fit the observed projections in different 
azimuths (Hanbury Brown et al. 1952). 


(e) Measurements of the Intensity Distributions across the Cygnus Source 
near the Major Axis 


The results for Cygnus, referred to above, did not appear to agree with 
the measurements made independently by Mills (1952) who concluded that 
the source exhibited a single central maximum with extensive limb 
darkening, the half-power width being approximately one minute of arc. 
The discrepancy could be explained if there was a dip in the correlation 
function in the neighbourhood of 1200), and a series of measurements was 
therefore made along the previous azimuth of 113° in order to investigate 
the intensity distribution across the source in a direction close to the 
major axis. 

A modification was made to the equipment in order to obtain the gradient 
of the correlation function, as the practical difficulties of obtaining small 
variations in aerial spacing were severe. The modification enabled the 
equipment to be operated on three alternative frequencies, 116-5, 125 
and 132-7 Me/s. The remote aerial and equipment were removed to a 
site where operation on these frequencies effectively provided baselines 
of 1450, 1558 and 1655, without the necessity of moving the aerials or 
equipment, and the assumption was made that the source distribution 
was not sharply dependent on frequency. 


of the Angular Diameter of Two Intense Radio Sources : II 69 


The correlation measured at these spacings was very low and the output 
of the cross multiplier due to the source was of the same order as the noise 
level of the equipment ; nevertheless it was possible to establish that the 
correlation function had a positive gradient. 

The equipment was then removed to a spacing of 2100 and the 
correlation coefficients were again obtained for three values of the effective 
baseline (Jennison and Das Gupta 1953). Further readings were later 
taken in the regions of the first maximum, the tail of the second maximum, 
and the locality of a possible third maximum. No third maximum was 
measurable although it must be remembered that the equipment was not 
suitable for the measurement of very small values of correlation and a 
small third maximum or tailing of the second maximum may exist. 

The readings obtained at the various spacings in position angle 113° 
are tabulated below, whilst the points are plotted for reference in fig. 1. 


Table 3. Cygnus Correlation/Aerial Spacing in Position Angle 113° 


Aerial spacing Correlation Aerial spacing Correlation 
in wavelengths coefficient (p?) in wavelengths coefficient (p?) 
610 0-460 +0-01 2230 0-290+0-01 
650 0-500 +0-05 2940 0-128 0-008 
900 0-200 +0-01 3140 0-085 +0-005 
1450 0-045 +0-005 3350 0-080 +0-01 
1558 0-065 + 0-005 4650 <0-025 
1655 | 0-120+0-01 5000 <0-028 
1952 0-245 +0-005 5300 <0-05 
2100 0-280 +0-005 
Rigel 
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(f) Further Measurements on Cassiopera 

Subsequent to the initial measurements, which served to confirm the 
identification, no attempt has been made specifically to investigate the 
structure of the Cassiopeia source. The post-detection interferometer is 
unsuitable for measurements of this type and its use is costly and cumber- 
some at short spacings. Readings were, however, taken concurrently with 
the Cygnus measurements at extended baselines in position angle 113° to 
enable a very approximate distribution to be derived. The values of 
correlation coefficients at these spacings are given in table 4, whilst the 
points are plotted in fig. 2 with the shorter baseline readings obtained by 
Smith (1952) at 210 Me/s in position angle 90°. The earlier establishment 
of a roughly circular distribution has been taken to justify the plotting of 
the points on a single graph, though this is only a temporary expedient to 
obtain the approximate distribution. In combining the readings in this 
manner it has been assumed that the source structure does not vary 
greatly between 116-5 Mc/s and 210 Mc/s and that the source is radially 
symmetrical. These assumptions are severe and further measurements 
are obviously required. 


Table 4. Cassiopeia Correlation/Aerial Spacing in Position Angle 113° 


Aerial spacing Correlation Aerial spacing Correlation 
in wavelengths coefficient (p?) in wavelengths coefficient (p?) 
605 0-103 0-004 2100 < 0-022 
650 0-092 +0-003 3140 <0-018 
900 0-050 +0-006 5000 <0-02 
1558 <0-02 
1.0) Fig. 2 
Oo. 
0.8 
SMITH 
0.7 ©  JENNISON & DAS Gupta. 
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Variation of correlation coefficients with baseline for the radio source Cassiopeia I. 
The plot is combined from the readings in position angle 113° given in 
table 4 and the measurements of Smith (1952) in position angle 90° 
Perfect symmetry has been assumed and no allowance has been made 
for the difference in azimuth. 
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§ 3. Discussion or REesuuts 
(a) Intensity Distribution of the Cygnus Source 

The correlation coefficients obtained for the Cygnus source suggest that 
the source consists of at least two prominent centres of emission. The 
curve could be generated by a simple symmetrical function or a complex 
odd function, and in the absence of phase information the discrimination 
between the alternatives cannot be made. A symmetrical two-body model 
of the source distribution which adequately fits the preliminary results is 
shown in fig. 3. A small amount of limb darkening has been given to the 
two sources in order to reduce the third maximum, which would otherwise 
have a value slightly in excess of the upper limit placed on it by the 
observations. A more complex solution is not proposed until this simple 
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(b) 


i i 31 istribution across the Cygnus 
Symmetrical two-body model of the intensity distribution 
: source. (a) Approximate overall distribution of radio intensity, the dotted 
contour represents the identified galactic collision. (b) The distribution 
of radio intensity in position angle 90°. 
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interpretation is disproved, for a more complicated model would in general 
have larger overall dimensions and would be even more difficult to reconcile 
with the astronomical evidence. 

In fig. 4 the measurements made by Smith (1952) and Mills (1952) are 
recorded on a single plot in which the above measurements have been 
converted into a position angle along the estimated major axis of the 
source. The full line in fig. 4 depicts the theoretical correlation function 
of fig. 3 along the major axis. It will be observed that the original 
discrepancy of the readings is largely explained if a model similar to that 
in fig. 3 is employed. 

Fig. 4 
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A combined plot showing the variation of correlation coefficients with baseline 
for the source Cygnus I. The tabulated readings for p? are given as 
ordinates whilst the abscissae have been foreshortened by projection on 
to the approximate position of the major axis. The readings obtained 
by Smith (1952) and Mills (1952) are shown for comparison. The full 
line shows the theoretical correlation function of the distribution shown 
in fig. 3. 

There is still some doubt concerning the minor axis of the distribution. 
Three measurements have been made in the north-south direction and 
these have yielded an effective strip diameter of 30”. As the width is 
relatively small a slight discrepancy between the line of measurement and 
the minor axis itself may lead to a considerable error of assessment, due 
to the fact that the interferometer measures the projected distribution 
in any axis. Thus the results could be equally explained if the minor 
axis 1s really 30” in position angle 180° or only a few seconds of are in 
a position angle a few degrees away from this azimuth. 


(b) Comparison of the Optical and Radio Frequency Distribution 
of the Cygnus Source 
The radio frequency distribution of the Cygnus source depicted in fig. 3 
is very much larger than the optical object and the position measurements 
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indicate that it must be presumed to straddle the position of the colliding 
galaxies. The dimensions are such that there can be little or no overlap 
between the regions of intense optical and radio emission and if the 
identification is correct, this must indicate that radio emission originates 
in the outermost and most tenuous parts of the galaxies, to the very limit 
where these might extend. The overall radio dimension in position angle 90° 
is 2'19", and at a distance of 310? parsecs this would indicate an 
actual width of 20400 parsecs which is just feasible for the overall 
diameter of the galaxies. The new distance-scale would double these 
figures, but the dimensions could still be consistent with the diameters of 
normal late type spirals. Ifany other model than the two-body distribution 
were chosen it would be necessary to consider a diameter considerably in 
excess of this figure. 

Ifthe visual identification and the interpretation of the radio distribution 
transform are correct, the most striking feature of the result is the negligible 
correlation, or possibly negative correlation, between the regions of optical 
and radio emission. This in turn provides an excellent test for any theories 
of the radio emission. Not only must the emission arise in the outermost 
parts but it must shun the bright regions in the centre. It may well be 
that later measurements will show finer details of the structure in the 
inner regions, but the presence of two very strong outer components 
appears to be indisputable. There is a faint possibility that the image of 
the source has been distorted by the interstellar medium en route to the 
earth, but the effect is unlikely to be serious. 


(c) Equivalent Black Body Temperature of the Cygnus Source 


The approximate evaluation of the Cygnus source distribution enables 
an estimate to be made of the equivalent surface temperature (7',) of the 
source from the known flux incident at the earth. The value of 7’, is 
given by 

SA? 
ad 2K fbdQ 


where S is the observed flux, 6 is the normalized brightness of the source 
and @ is an associated element of solid angle. If the first order distribution 
of the Cygnus source is taken as consisting of two rectangles 50” x 30", 
no limb darkening being assumed, the apparent temperature at 125 Mc/s 
is 71,3108 deg. x. This figure is obviously higher than the thermal 
temperature which would eventually result from the collision. In the 
initial collision the electron energies are small and, in the absence of 
radiation, the proton energies may be taken as approximately representing 
twice the final particle energy on equipartition. Thus for a collision at 
1500 km/sec, 7,~5x10’ deg. x. If the mean density is taken as 
1 atom/cm? this corresponds to an absorption coefficient seth parsec™* 
and the galaxies viewed face-on should therefore appear optically thin. 
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If thie same source distribution still holds at very much lower frequencies 
the discrepancy is even greater. An apparent temperature ~10" at 
30 Mc/s can hardly be reconciled with a simple thermal process at the 
source. ' 


(d) Intensity Distribution of the Cassiopera Source 
The most prominent features of the curve in fig. 2 are the long tail and 
the absence of pronounced subsidiary maxima. A curve of this form 
would result from the transform of a tapered distribution and if a 
symmetrical distribution is assumed the source may be considered as 
having severe limb darkening whilst the long tail could be generated by 
a smaller intense region situated near the centre of the general distribution. 
This model is shown in fig. 5 but the interpretation must be treated with 
considerable uncertainty in view of the small number of useful observations 
and ‘the!assumption of symmetry. There may well be marked deviations 
from the,curve in fig. 2 when more readings have been established, and the 
distribution in fig. 5 will have to be modified accordingly. Ifthe true curve 
is complex it is also probable that the phase of the transform will indicate 
that a symmetrical distribution is not justified. 


Fig. 5 


Symmetrical model of the intensity distribution across the Cassiopeia source. 
This model is given only as a preliminary result and may require extensive 
modification when investigated with a more powerful technique. 


(e) Comparison of the Optical and Radio Frequency Distribution 
of the Cassiopeia Source 
Some justification for the symmetrical model of the Cassiopeia dis- 
tribution shown in fig. 5 is provided by the position of the source, The 
accurate radio position, as measured by Smith (1951) coincides with the 
centre of the region as a whole, whilst the photographic centre of gravity 
is considerably removed towards the arch on the northerly limb. The 
photographic distribution differs in the blue and red, and, as in the blue 
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only the arch and a few diffuse filaments are visible, whilst in the red the 
arch becomes less prominent and more filaments appear over the whole 
area, it 1s apparent that a detailed comparison between the optical and 
radio frequency distributions must take into account the passband of the 
photographic filter, whilst it is also possible that the radio frequency 
distribution may differ at the extremes of the radio spectrum. 

Though the accuracy of the present measurement is far too low to permit 
a detailed analysis, it is apparent that, as with the Cygnus source, the 
correlation between regions of optical and radio emission is poor. It is 
evidently of great importance to investigate the radio distribution more 
thoroughly, and though the resolution required to discern detail of the 
size of the optical filaments may prove extremely difficult to obtain, the 
whole assemblage of filaments in the northerly limb may be more simply 
resolved from the main body of the nebula. 


(f) Equivalent Black Body Temperature of the Cassiopeia Source 
The approximate intensity distribution of the Cassiopeia source may be 
used to derive the equivalent temperature of the source. The resulting 
apparent temperature: at 125 Mc/s is ~4 10" deg. kK over the main 
body of the structure. 
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ABSTRACT 


The Kittel-Devonshire phenomenological treatment of antiferro- 
electricity is extended to the case of polarizations in three dimensions 
in a pseudo-cubic material. There are ten possible solutions of the 
_ proposed free energy equation, corresponding to non-polar, ferroelectric, 
antiferroelectric and ferrielectric states, and expressions are derived for 
the polarization and dielectric stiffness coefficients in these forms. 

Using the experimental data for NaNbO,, values of the constants im 
the free energy expression are derived which give satisfactory agreement 
with the observed phase changes, optical anisotropies and dielectric 
constants in this material. It is also shown that the major features of 
the phase diagram for (NaK)NbO, solid solutions are reproduced from 
the NaNbO, free energy expression if 7',, the antiferroelectric ‘ curie ’ 
temperature decreases linearly with increasing KNbO, concentration. 
These results are discussed in terms of the spatial relationships in the 
NaNbO, and KNbO, structures, and compared to the properties of 
the corresponding alkali tantalates. 


§ 1. INTRODUCTION 


THE thermodynamics of antiferroelectricity was first studied by Kittel 
(1951). In this treatment, the dielectric is assumed to contain two 
equivalent, interpenetrating, polarizable sub-lattices. Designating the 
sub-lattice polarizations P,, P,, the free energy is expressed in the form 


AG=f(P2+P)+9P PtP s+PAt+.... . . (1.1) 
There are two possible solutions of (1.1), P,=P,, a ferroelectric state, 
and P,=—P,, an antiferroelectric state. The stable solution is that 


which gives the least minimum value of (1.1). 
Devonshire (1954) shows that it is more convenient to transform (1.1) 
by putting in new variables 


Pat PtP or pee aa yy meta ’ 
EER aR eee te) Del eed ey cle) 
The free energy expression then takes the form 
AG=3x'P?+ hE'PA+.. .4yp2+fépt+...+hrpePet.... . (1.3) 


* Communicated by the Author. 
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The conditions for stability are 


CAC OL == Onin en es 1. 2 1i4) 
CACC UO. le er eee) 3. (15) 

and the dielectric stiffness is given by 
Y=O-AG OL, er veneer ei .G) 


It has been suggested (Kittel 1951, Shirane et al. 1951) that IMeyAcO)s. 
and certain related compounds are antiferroelectric. The thermodynamic 
method gives a plausible explanation of the dielectric anomalies at the 
curie temperature in these substances, but its usefulness is limited by 
lack of experimental data which prevents any detailed evaluation of the 
coefficients in (1.3). 

Recently, more complete experimental data have been published for 
NaNbO;, and mixed (KNa)NbO,, crystals, and it would appear that a 
more detailed thermodynamic analysis might profitably be attempted. 


§ 2. DEVONSHIRE’S EQuaTION FOR A PsEUDO-CuBIC MATERIAL 


In (1.3) both P and p are vectors, and in general they need not lie in 
the same direction. For a pseudo-cubic material (1.3) may be expanded 
in rectangular coordinates, and takes a reasonably simple form. 

Assuming again that the material contains two interpenetrating 
polarizable sub-lattices with polarizations P,, P,, these will have 
PMiIpOnciiselo nn! 7 and PoP. P.. 1m the directions. of the 
cube axes x, y, 2. 


Putting 
ete a Pa=Poz—Py, (4) | 
Py=P,,+Pr, (6) ls Oe (2.1) 
iy een) Sete cad 4 


the free energy may be expressed in the form 
AG=3x(PP+PP+P2)+ 4b (Pott Pitt PA) 

pereieet  lyle te ) ell, +e, +e) 

+4Xa(Pa? + PP t+De)+ Fo (Po + Py + Pe) 

+4£12( D2" t+ Py Pe +De Pa) + obra + Pi TPs) 

+ Sire Pe Py? + Pe POPs Pa + Py Ps + Pe Pe t Pe Py’) 

+ Bip pP + hy (Pert Py tPA) tah eps + Py py + P.'P.") 

440(P2p2+Pp2+P 2p 2+Pep2+P2p2+P2p,y). . - (2.2) 

The ‘ ferroelectric’ part of the expression contains all possible terms 

up to the fourth order, and the first sixth order term. The " antiferro- 
electric ’ part contains all terms up to the sixth order, and the first eighth 
order term. Only first order interactions between ferroelectric and 
antiferroelectric are considered. 
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From (1.4) the conditions for stability in zero external field will be 
P,=0 or 0= yy Eqy'P 2+ £49'(PP+P2)+ OP gt + Ape +O(P,2+P2), (a) 
P,=0 or O= yt E41'Py?+ E12 (Pa2 + P2)+OP A+ Ap + (ps +P), (6) 
P,=0 or O= yy +Ey1'P2 +810 (Pet PY)+oP A+ Ape + Ops + Py), (C) | 
and from (1.5) pee 

Pr=9 or O= xg tF1rP 2? + 10(Ps? +P.) + S111P 2 | 
+ bbira{4(De Py + PsP?) +2(pyt pe)s | 
+a Py Pet IDe Al gt lg ole Ae) | 
y= or O= Xqt yy? + E12 Pa + Pl) tory | 
+ $b 112{4( Dy D+ Dy Pa") +2(p2'+pe')} (2.4) 
+ ees P oP + py +AP +P +PF), (0) | 
p,=9 or O= x4 + EPP + $12( Pe + Py) toi Ps 
+ $Sarot4(De Ds + PPD") +2(Dx + Py’) } | 
yes Pe Oy ees ee eee J 
From the symmetry also, 
if 2,960, then) p7==0), =) = a ees 
and if », 340, then] P2070 =. ee 


where n=, y, 2. 
There are then ten possible sets of solutions for (2.2). 


Pe ive tie P2=Py=P,—9, (a) 
O= xt (En )P2+EP A, P,=P,=0, 
Pr=Py=P-=9, (6) 
Omer (E11 4 £19 hg Rees P,=P,, P,=0, 
Pr=Py=P2=9,  (C) 
Oe Nets (Oxy te eae eget Ces P, =P, =P,, 
Pr=Py=p),—9, (d) 
fay ead eA 
O=XgtF11Pe2teiiPs t+npe's Py=p-=9, (e) 
Aes Wai wage IN), 


O= Kat (S11t £12) Part (S131 +5112) Po? +728, 
Pre=Py p,=9, (f) 
Pf P90, 
O= Xat (S1rt 2819) Pa7 + (S411 + 201194 S303) 2? + p28, 
Pr=Py=P» (9) 
O=xX+611'P2+0'P,'+6p,?, igh Ve ANY 
O=xX.+'110," +b, END sy 640P.2, *, Pz=P,=0, (A) 
O=xX+ (E11 + £10 )P.2+0'P 4+ Op,?, P—P,, P=, 
O=Xot F11P2 + 01110-4-+-np,S+ 20P,,2, Pr=P,y=9, (2) 
O=yxy- +617 P74+0P A+ 20D"; hes P i= 0, 
O=Xot (Sri tE19)p,7 (Gs yoeC ies wo UTE ean cee 
P:=Py, Pr=9. (9) 
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The stable solution at any temperature will be that which corresponds 
to the least minimum in the free energy (2.2). 

From (1.6) the dielectric stiffness coefficients are given by 
Kaw — Nyt 8841 Pt Ess (PP + PZ)+50P, 4 Ap 2" O(p,? (oe) (a) 
Xuv =X 3611 Py? bya (PP +P?) +5C'P,*+-Ap,2+-6(p,2+-p.2), (b) (2.8) 
Kea = Xo 8611 PP + bie (PP + Py) +50 P44 Ap? + O(p,2+p,2). (6) | 


§ 3. APPLICATION OF THE THERMODYNAMIC TREATMENT TO N aNbO, 


To make use of this treatment in NaNbOsg, it is necessary to determine 
the constants in both ferroelectric and antiferroelectric parts of the free 
energy expression (2.2). The situation is thus more complex than in 
BaTiO, and it is necessary to make more direct use of the experimental 
data than in Devonshire’s treatment (Devonshire 1949). 


3.1. Hxperumental Data 
The necessary data are taken from papers by Cross and Nicholson (1955) 
for NaNbOg,, Shirane, Newnham and Pepinsky (1954) for NaNbO, and 
mixed (NaK)NbO,, and Shirane, Danner, Pavlovic and Pepinsky (1954) 
for KNbO,. In further references in the text, these papers are abbreviated 
C,, 8S, and S, respectively. 

Qualitatively, the following sequence of properties with decreasing 

temperature is suggested : 

(a) T>638°C Cubic probably paraelectric. 

(6) 638°c>T7>500°c Tetragonal, almost cubic, the distortion being 
too small to detect by x-rays, probably para- 
electric. 

(c) 500°C >T7'>350°C Orthorhombic, almost tetragonal, again the 
distortion from the higher symmetry may only 
be detected optically, probably antiferro- 


electric. 
(d) 350°c>7>—135°c Orthorhombic (rotated axes) antiferroelectric. 
(e) —135°c>T Probably monoclinic, probably ferrielectric 


with antiferroelectric displacements in the 
a.c. plane, and a ferroelectric polarization 
along (b). 

There appears to be some doubt as to the properties near 500°C. Cy 
reports minor structure changes at 470°c and 518°c, 5, a minor change 
at 480°c. Evidence for the low temperature polarized state below —135°c 
comes only from C, ; 8, reports no symmetry changes at low temperature 
but does not appear to have investigated this region in detail. 

The general data suggest that the sequence of solutions of (2.2) with 
decreasing temperature is (2.7(a), (e), (f), (9): and these solutions will be 
used to evaluate the constants in the free energy expression. 
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3.2. Evaluation of Constants in Free Energy Equation 


Following Devonshire 1949, it will be assumed that x, is a linearly 
decreasing function of 7', of the form 


Xs=Xr0(T—T,) a ee 
and also that y,, has the similar form 
Xa=Xaoll’—T,) ey o> tie) Rome eaaes (3.2) 


where y,, 7';. X49 and 7, are constant. 

In an antiferroelectric, the macroscopic polarization is always zero, so 
that it is not possible to measure p directly. It is probable, however, 
that p will be related to the structure distortion, and to the optical 
anisotropy, and it will be assumed that this relation takes a similar simple 
form to that found for the ferroelectric polarization in BaTiOg, i.e. 


v2=p(z—1) SKA Oe ee ee 


where p and K are constant. 

In the orthorhombic form there are two measurable anisotropies. For 
this work Au, (C,) will be used, as it is the difference between two principle 
indices. 

From (3.2) and (3.3) solutions (2.7(e), (f)), of (2.2) may be rewritten in 
terms of Mu and du, 


0=(T—T,)+AAu+C Ap?+ F Ap, ee «| ee ee een 
0=(P—T,)+(AFB)Auy+(CHD) Au e+ PAu? . . (3.5) 
and from (2.2) the free energy in the two states will be 
AG=K yao {3(T—T)4u+}4A Au?+2C Ape +ifdpt} 2. 2. . . . (8.6) 
AG=K Yo {(T—Lq) 44+ MA+B) An P+HC+D) An 3+}F 4p}. (3.7) 
where 


A=KE41/xXa0; B=KE12/X a0; C= KC 11/Xa0, D=K%E112/Xa0> #=K*n/xX,0- 

To evaluate the constants A, B, C, D, F the following procedure is 
adopted : 

(a) By the method of least squares, values of 7',, A+B, C+D, and F 
are obtained which give a * best fit ’ between (3.5) and the experimental 
temperature variation of 4u,, in the orthorhombic state (C,). 

(b) At the transition temperature from orthorhombic to tetragonal 
structures, the free energies of the two states must be equal, i.e. 

2(2— Ty) Aut fA Ap? + C Ape + fF Apt 

= (2, —Tq) Au, +}(A+B) Apu +3(C+D) 434 4h Ay! . (3.8) 
substituting the experimental values of 7,, du, du, and the constants 
previously determined in (a), (3.4) and (3.8) may be solved for A and C. 

The least squares method (a) gives a value of 7.,=538°c ; this is close 
to the observed minor structure change at 520°c, and assuming this to be 
the antiferroelectric * curie’ temperature recalculation with 7',—520°c 
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‘gives the values A+ B=1-322 x 104, C+ D=—2-526 x 10°, F=2-326 x 106 
and the fit with the experimental 4, shown in fig. 1. 
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The experimental value of Jy at the transition temperature substituted 


in (6) gives values of A and C, which in (3.4) produce a temperature 
variation of Jy dissimilar to the observed variation (fig. 2, curve (a)). 
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Using a value of 4.—0-034 at the transition, a more similar temperature 
variation is obtained (fig. 2, curve (b)). The difference in scale between 
observed and calculated values may be due to the oversimplification of 
using a single value for K in both orthorhombic and tetragonal regions. 

The values of A and C for 4=0-034 are 

A=0-612>¢10% C=—0-989 x 10°. 

Using the constants enumerated above, the free energy in units of Kyo 

may be evaluated from (3.6) and (3.7), and is shown in fig. 3 (curves (@), (5)). 


Fig. 3 
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Free energy in NaNbO, as a function of temperature. 


3.3. The Determination of the Ferroelectric Constants, 
and the Cross Constants 


The ferroelectric curie temperature 7’,, x; and the cross constants ) 
and 6, may be determined from the dielectric constant data. 
For solutions (2.7(a), (e), (f)) eqns. (2.8) may be rewritten 


T'>520°c Xax—Xyy= Xoo X02 — T's) (3. (3.9) 
520°C>T'>320°o - xy ,aeeng( let) haan eee WReray 

Xvv =Xee=Xr0(T T;)+ BAu Ghac Ly) 
350°c >7' > — 209°C Xea=Xvy=Xp(2—T,y)+(a+B)4u, . . (a)) 

Kez =Xa(T—Ty) +28 Ap, t @pew 


where poe 
eek and peeks. ov. (3:12) 
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The values of the constants yj, 7',, « and B giving the best fit between 
calculated and observed dielectric constant are Mea 0; CLOSE == 3 10°C, 


%=1-37, B=0-265 and the fit obtained is shown in fig. 4. 


Fig. 4 
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To determine the remaining ferroelectric constants, it is necessary to 
consider the low temperature ferrielectric state. Experimental data are 
meagre, but the sudden onset of hysteresis and the small change in 
extrapolated polarization with temperature (C;) would suggest that at 
the ferroelectric curie temperature there would be a first order change. 
With this assumption we may make use of Devonshire’s treatment (1949) 
to simplify the ferroelectric part of (2.7(j)). Substituting also from (3.3) 
for p, this becomes 

0=((T—T,) (14-1) + 28 4p y/xXp0(L1—T's)) 42 +32 . (3.13) 
where 
7’, is the ferroelectric transition temperature, 
Ay, a measure of the optical anisotropy in the ferrielectric, 
P is the ferroelectric polarization at 7',, 


P\2 
z= (—* Re ere es EE | 
By a 


The * antiferroelectric ’ part of (2.7(7)) becomes 
0=(T—T,)+Bwel {xio(L1—T)}4+ (A+B) py + (C+D) Ap? +F Ap? 
oath (3.15) 
where W= VX polly) pe Exe) 
and the free energy in the ferrielectric state is given by 
AG=Kyqn{(GP—2-+-42(T —T,) (Ps —T,) 0+ (T—Tg) Aye + HA+B) Mpeg? 
44(CHD) Ap yP+ EF Ap yt + BApytee|(x(P1—Ly))}- + (8.17) 
At the transition temperature 7',, antiferroelectric and ferrielectric states 
must have similar free energy so that from (3.17) and (3.7) 
(32° — 27+ 32(T'p--T) (TT) wt (Tp9—T a) Ati 
+3(A+B) Mpg? +3(C+D) Ap ye + EF Apa +B Ap js w2/xX4o(T 1—7's) 
=(Tp—Tq)Apy+3(At+B) Ap 2+3(C+D) Au 3+ fF Aes (3.18) 
With the additional condition that no other ferrielectric or ferroelectric 
states have lower free energy than the antiferroelectric structures observed 
at higher temperature, a reasonable fit with the observed transition 
temperature i obtained with 
w=3-65, 1,=415°. 
The free energy in the ferrielectric phase, using these values is shown in 
fig. 3 (curve (c)). The alternative ferroelectric tetragonal phase (fig. 3, 


curve (d)) deduced using these values of wand 7’, has only slightly higher 
free energy than the antiferroelectric states at 335°c (fig. 3, curves (a) (b)). 


3.4. Dielectric Constant at Low Temperature 
For the ferrielectric solution (2.7(j)), the dielectric stiffness eqns. (2.8) 


become 
Xxx Xve=Xp0(1'—T,)+(A+6) p72 +10 jet eG ONS) 
Xao= Xy0(2'— 7 ,)+364)'P,?+50'P 44 20p,?, . . (3.20) 
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Using eqn. (3.2) and Devonshire’s simplification for the ferroelectric 
terms, (3.19) and (3.20) become 

Xan Xyv=Xso( LT —T) + (x +B) Aa t+ 4x4o(2'1—T az, - + (3.21) 

Xe=Xp0(01—T) {((T —T,) (0, —T ,)— 1224 1527} + 28 Aug. (3.22) 


Substituting the values of z, Au, and 7' calculated in § 3.3, eat erp ee 
may be evaluated, and thus the corresponding dielectric constant values 
deduced. These are shown in fig. 4(b). 


3.5. Hlectrically Forced Transitions 

The experimental double hysteresis curves at room temperature (C,) 
show that there is a sudden increase in the total polarization at a field 
100 kv/em, applied in a direction perpendicular to the orthorhombic 
c-axis. 

From (2.2) there is a solution (2.7(/)), which is a spontaneously polarized 
ferrielectric state of only slightly higher free energy than the stable anti- 
ferroelectric orthorhombic state at zero field. At high field strength this 
may become the state of least energy. 

For a field #, applied axially along x, conditions for stability of (2.7(h)) 
become 

Pee oro 44 lea O Lacey), 8s =. + » ~ (3.28) 
UN oii y ity Py Ol : : 9 (3.24) 

Using (3.3) and Devonshire’s simplification for the ferroelectric terms, 

(3.23) and (3.24) may be rewritten 


E,| {xP (1 1—T) =? {(T—-T 7) (11-7) +B Ap / X(T 1-7) — 42+ 827} 
: (3.25) 
0=(P—Ta) 4 Bwe|(xp(Pr—T)) +A Au +-CAwe+F Aue — (3.26) 
and from (2.2) the free energy is given by 
AG= (42-24 4e(T —T,)|(T,—T,) w+ MT —Ta) Apt $4 Ap? 
+40 AP + 4F Aut + $8 Auwe|(x5(L 1 —T)) +B 2*?0/P 1X 0o(L1—1'}). 

(3.27) 
The dielectric constant in the antiferroelectric state is low for this field 


orientation, so that for an approximate solution we may neglect the energy 
contribution from the field in the antiferroelectric state, and write 


0=(T—Ta)+(A+B)Auy+(C+D)Aue+F Ans,  . . (3.28) 
AG=(T—Ta) Any +H A+B) Ap 2+HC+D) Mn 2+4F Au). (3.29) 
At the forced transition 
—22+-12(7—T,)|(1,—T,;) }wt oT T ,)Au,+4A Ap? +46 Ape EF Apt 
AB Apwel(x4o(L1—Ty)) + E ee] xpoP (71-7 )=(F—-T a) Aen 
AA(AERB) Ag Pt HO+D) Ape teh Ag 9.2. + ee 13.80 
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Eliminating #, and z, (3.30) may be solved numerically for du at any | 
fixed temperature, and substituting back in (3.25) and (3.26) z and H/P, 
evaluated. 

Using P,=10,.c/em? the calculated temperature variation of the field 
required to force the transition is shown in fig. 5. 
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$4. THE Properties oF MrxEp (KNa)NbO; CrysTALs 


KNbO, has been shown to have properties closely analogous to those 
of BaTiO, (S,). There are first order structure changes at 425°C, 225°c 
and —10°c, and the observed symmetries suggest that the stable states 
correspond to solutions (2.7(a), (6), (c), (d)) of (2.2). 

The properties appear to be adequately described by a simple Devonshire 
expression with values of the constants 

T ,=430°c, 1,=320°C, x 4o=2°62X 10-5, a=0-365. 

In NaNbOs, 7’, and 7', have similar values, respectively, to those for 
KNbO,. Assuming that « is also similar for the two materials, the free 
energy in the unstable* orthorhombic and rhombohedral states may be 
calculated (fig. 3, curves (e), (f)). 

It has been shown (§,) that in mixed (NaK)NbO, crystals the anti- 
ferroelectric ‘curie’ temperature is linearly depressed with increasing 
KNbO;. If in the free energy expression (2.2) only 7’, changes, the free 
energy of mixed crystals may be directly deduced from fig. 3 by moving 
the antiferroelectric curve to lower temperature without changing its 
shape. 


a Ee 

* The ferroelectric states are described as unstable in pure NaNbO, as though 
the solutions (2.7 (b), (c), (d)) correspond to minima in the free energy, they are 
not at any temperature the least minimum. . 
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Using the experimentally determined value for the rate of change of 
curie temperature with concentration, i.e. 
T= 10°c/1% added K in (NaK)NbO,, 


the free energy for a number of concentrations of KNbO, is shown in 
fig. 6. From these and other similar curves, the phase diagram for the 
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system NaNbO, to (Nag.gKy.;)NbO; has been deduced (fig. 7). The 
behaviour of the antiferroelectric to ferrielectric phase transition is not 
easily determined graphically, and the values shown in fig. 7 have been 
calculated from (3.13), (3.14), (3.15), (3.16) with suitable values of 7',. 
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§ 5. Discussion 


5.1. Comparison with Experimental Values 


Over the orthorhombic temperature range (— 200°C to +350°C) there 
is reasonable agreement between the observed and calculated anisotropy 
values, and as has already been indicated, the ‘ best fit ’ between (3.5) 
and the experimental curve gives 7’, close to one of the observed transitions. 
To obtain this agreement, it is necessary to include the term in Ap? in (3.5) 
and thus to include terms up to p® in the free energy expression. It should 
be noted also, that though the constant ¢,; cannot be evaluated from 
the available data, K2Z,53/x,9 must be greater than 2-035 x 10° to avoid 
the occurrence of a rhombohedral antiferroelectric phase above —200°c. 
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To obtain agreement with the optical anisotropy in the ‘ tetragonal ’ 
region, it is necessary to assume different values of K, the constant in (3.3), 
in orthorhombic and tetragonal regions. From (3.3) however, there 
should also be a linear relation between p? and the elastic anisotropy, and 
consideration of the electrostrictive constants suggests that the constant p 
in (3.5) should be the same in both phases. Plotting c/a—1 against 7’ (S,) 
in figs. 1, 2 does show increased anisotropy in the tetragonal region, but 
the agreement is still poor. 


of Ferroelectricity and Antiferroelectricity 89 


The free energy expression (2.2) assumes that in the unpolarized un- 
stressed condition NaNbO, has cubic symmetry. Optical measurements 
(C,) however show that above 520°c the symmetry is not cubic, and below 
520°c the dielectric constant values are not consistent lis tetragonal 
symmetry. Clearly (2.2) is only an approximation to the true free energy 
expression which must have lower symmetry, but in view of the complexity 
of the expression no attempt will be made here to improve the agreement. 
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Phase diagram for mixed (NaK)NbOs. 


In the thermodynamic treatment it has been assumed that the low 
temperature polar phase is ferrielectric. The general agreement with 
observed properties would suggest that this assumption is essentially 
correct. A critical test could, however, be applied, as in the low tempera- 
ture form, x-ray examination should reveal the super-lattice associated 
with the antiferroelectric polarizations. 

A striking feature of the mixed (NaK)NbO, crystals is the rapidity 
with which the properties change from those of pure NaNbO,; with in- 
creasing KNbO, concentration, so that in solid solutions with less than 
90% NaNbO, the properties are very similar to those of pure KNbO,. 
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In the foregoing analysis this rapid change arises through the presence 
of latent ferroelectric forms in pure NaNbOg, similar to the stable phases 
of KNbO, but of slightly higher free energy than the antiferroelectric. 
The general features of the phase diagram for mixed crystals are well 
reproduced from the NaNbO, free energy expression changing only the 
parameters 7), and xzo. 

Below room temperature, the data for mixed crystals are not yet 
complete, but there is no evidence for the rhombohedral ferroelectric phase 
which appears in the analysis for KNbO, concentrations greater than 6%. 
The treatment however also predicts that the ferrielectric—antiferroelectric 
phase transition temperature should rise with increasing KNbO, concen- 
tration in the mixed crystals. It is possible that the chosen values of w 
and 7', put this phase change in pure NaNbO, at too low a temperature 
(the exact transition temperature is difficult to confirm experimentally 
because of a wide thermal hysteresis) and that values could be chosen 
which would make the ferrielectric more stable than the rhombohedral 
ferroelectric phase for intermediate compositions. Further experimental 
evidence is necessary to decide between these alternatives. 


5.2. Correlation with Crystal Structure 


The phenomenological method gives no direct information as to the 
origin of the observed properties, and the necessary complementary 
statistical mechanical treatment will not be attempted here. Certain 
features of the results with the mixed crystals do however suggest a 
correlation with the spatial relations in the structure. 

In pure KNbOg, the radii of the ions are such that the Goldschmidt 
tolerance factor ¢ (Goldschmidt 1927) for the perovskite structure is very 
close to 1, i.e. the K ion is only slightly too small to fulfil the requirements 
for close packing. The structure may thus be regarded as a regular array 
of corner linked NbOg octahedra with the K ions almost filling the inter- 
stices. Substituting the smaller Na ion in some of the K sites in the 
(Nak)NbO, solid solutions will clearly have little effect on the already 
close-packed octahedra. ‘The very small change in the ferroelectric 
properties with increasing NaNbO, would thus suggest that in the solid 
solutions the ferroelectric phenomena are essentially determined within 
the NbO, array. 

At high NaNbO, concentrations, the regular NbO, array may become 
unstable, “ crumpling ’ to reduce the excessive space available to the Na 
ions. A possible explanation of the high temperature structure changes 
in NaNbO, may then be given as follows : 


(a) At 638°C the regular NbO, array ‘ crumples ’ to reduce the excessive 
space about the Na ions. The symmetry is reduced to tetragonal with 
cla<1, 


(b) At 520°C ordered displacements occur in the new NbOg array, 


alternate sheets of Nb ions perpendicular to the tetrad axis being displaced 
im opposite sense. 


of Ferroelectricity and Antiferroelectricity 91 


(c) At 470°c the character of the ordering changes by a doubling of the 
repeat interval, i.e. the Nb ions in successive sheets perpendicular to ¢ 
have polarizations -- ++ ++ ete. 


With this sequence of changes, the displacements of ions occurring in 
the room temperature phase (Vousden 1951) may be accomplished at 320°c 
without a major reconstruction of the existing system. 


5.3. Properties of the Alkali Tantalates 

Some support for the dependence on spatial relationships proposed in 
§ 5.2 may be obtained from a comparison of the properties of the niobates 
with the corresponding tantalates NaTaO, and KTaO. Tantalum and 
niobium have very similar ionic sizes, so that properties which are essentially 
“space dependent ’ should be common to the two structures. 

That the onset of ferroelectricity is not primarily determined by the 
space available to the B ion in these ABO, compounds is clearly demon- 
strated by the potassium salts. In KNbO, the curie temperature is at 
425°c, whilst in KTaO, in spite of a slightly smaller Ta ion, spontaneous 
polarization does not occur above 14°k (Hulm and Matthias 1950). 

In sodium niobate, it is suggested in § 5.2 that the collapse of the structure 
about the Na ions which initiates the antiferroelectricity is primarily due 
to spatial considerations. We have recently investigated optically the 
properties of single crystals of NaTaO, at high temperature (Cross 1955) 
and do find a similar transition from cubic to tetragonal symmetry at 630°C. 
There are also further phase changes to an orthorhombic (pseudo-tetragonal) 
state at 580°c, and the room temperature orthorhombic (rotated axis) form 
at 480°c. (This latter transition would appear to be the change previously 
reported as the curie temperature (Matthias 1949).) There is no evidence 
of any change corresponding to that occurring at 470°C in NaNbO, and it 
is interesting to note that a recent determination of the NaTaO, structure 
at room temperature (Kay Miles 1954) shows a doubling of the b spacing 
of the primitive unit cell, as compared with the quadrupling observed in 
NaNbOs,. 

Unfortunately there are no data regarding the dielectric constant in 
NaTaO, and no information appears yet to have been published on mixed 
(KNa)TaO, solid solutions. 

List of Symbols 
4G The Helmholtz free energy of a free crystal per unit volume. 
(Taking the energy zero to be that of the unstressed 
unpolarized crystal.) 


dente The sub-lattice polarizations. 

P The total polarization (P,+P,). 

p A measure of the antiferroelectric polarizations (P,—P,). 
1 ced eit ‘The components of total polarization and antiferroelectric 
ee Day Ds Js polarizations in the axial directions. 


» tare The coefficients of the ferroelectric terms in the free energy 
Sas 6. "expression. 
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Xa» €11> £12 S111 The coefficients of the antiferroelectric terms in the free 
C112» S123) 7 energy expression. 


A, 6 The coefficients of the ferroelectric, antiferroelectric cross 
terms in the free energy expression. 

Xr0 The temperature invariant factor in y,. 

Xa0 The temperature invariant factor in x,. 

ee The ferroelectric extrapolated curie temperature. 

f he The antiferroelectric curie temperature. 

Y The dielectric stiffness, i.e. the reciprocal elastic suscepti- 
bility. 

Oey yuk Nae The axial components of the dielectric stiffness. 

An The birefringence of tetragonal NaNbO3. 

Ap, The birefringence B—« of orthorhombic NaNbOs3. 

Aum A measure of the optical anisotropy in the low temperature 
form of NaNbOs. 

p The constant of proportionality between elastic aniso- 
tropy (c/a—1) and the square of the antiferroelectric 
polarization. 

Kk The constant of proportionality between the optical 
anisotropy and the square of the antiferroelectric 
polarization. 

Te The temperature at which the unstable ferroelectric and 
paraelectric states have the same free energy. 

Le The spontaneous polarization in the ferroelectric state 
at T,. 
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VIII. An Unusual Double V-Event at Sea Level 


By J. L. Luoyp and A. W. WotrenpaLe 
The Physical Laboratories, The University, Manchester* 


[Received September 27, 1955] 


A MULTIPLATE cloud chamber has been operated at sea level under the 
Manchester Spectrograph in an experiment designed primarily to study 
the scattering of y-mesons (Lloyd and Wolfendale 1955 a) and among the 
photographs taken is the one reproduced in Plate 1 and the figure. 
The spectrograph gives no information on the momentum of the particle 
responsible for the interaction since the primary particle is part of a 
shower. 

The event has been reprojected, through the actual lenses used, 
together with other photographs taken before and after in order to find the 
extent of the illuminated region, condensation and illumination conditions 
and noise-level scattering. The photograph probably indicates an 
interaction by the particle marked (1) in one of the 1 cm lead plates, in 
which is produced a number of charged particles, (2)-(5), together with 
two neutral particles emitted in opposite directions and decaying into 
particles a, and a,, 6, and b, respectively. The knock-on electron 
produced by (1) in the top plate shows that this particle is very probably 
moving downwards. 

The primary particle has the ionization of a fast singly-charged particle 
(Jat) and an r.m.s. projected angle of scattering of 0-4° per plate ; tracks 
in this region have an r.m.s. noise-level scattering of 0-2°. Assuming 
particle (1) to be a proton the estimated momentum is 2-1+7? eev/c 
(Lloyd and Wolfendale 1955 b). The most probable value is quoted 
and the errors are the limits for a 68-27°% probability of the true value 
lying between them. If (1) is a7-meson the momentum will be approxi- 
mately the same. The data of Camerini ¢¢ al. (1952) on star production 
indicate a momentum of about 4 Gev/c, in rough agreement. 

Particles (2), (3) and (5) pass out of the illuminated region with 
ionizations less than 3X Lpjat- Particle (4) has ionization (1-5-4) x Iyiat- 
From the scattering in two plates the value of pf for particle (5) is found 
‘to be p8=1407% mev/c. From considerations of range particle (5) is 
probably a 7-meson although a p-meson is not excluded. 

The origin of the interaction is well defined by tracks (1), (2), (3) 
and (5). 

Two neutral V-particles, designated a and b are associated with the 
interaction. Particle a decays into two charged secondaries a, and ay, 
8-4 em from the origin of the interaction. The spatial angles ¢, and 4, 
between the trajectories of a and a,, a and a, together with the estimated 


* Communicated by the Authors. 
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ionizations are given in table 1. Evidence that the decay event is 

directly associated with the interaction is provided by the fact that the 

lack of coplanarity (the angle 5 between the plane of the secondaries 

and the origin of the interaction) is small, 6 being 1-5+-2°. Particle a, 

penetrates two plates before reaching the edge of the illuminated region. 

Scattering measurements indicate that the momentum is given by 
pB=800 *$8° Mev/c. 


Edge of 
illuminated 
region 


Side View 


Particle a, travels obliquely towards the rear and the track would 
not be seen to emerge even if the particle had penetrated the plate. 

The upper V-event consists of two decay products b, and by, the lack 
of coplanarity again being small: 6=1+-6°. The track of particle b, 
would be visible above the plate if the particle had penetrated (unless it 
had suffered a large deflection towards the front of the chamber), so 5, 
has a range less than 1-7 cm Pb. 

The V-events have been analysed (following Podolanski and Armen- 
teros 1954) on the assumption that they represent the dacays of 6°- or 


95 


$0LFS SI USual Avoap uvout 07 YySuet yyed paadesqo Jo orea IuL 


Unusual Double V-Event at Sea Level 


a a eee a ae a a a 
rae S*P6 (4) @-EFS-06 C-3-G % 
eS 6 oF 006 (d) 8F 8s 6-F "@ 

(o]B9) (-oT@o) 9/Aew eek, 
Lael ae UL WH {USULOP, ‘ ne 
; (Sap) UOT} eZIUOT apolaed 
, oF p SATYLIOI JO A.epuooeg 
9/AOW BF LLToV O}BUITYSO TRUST A 
aq 04 ties: Areutlig : ; 


9 TU9AH oA | SG PIVEL 


‘ouloyos ABoop Yowe JOF AZIUN Jnoqe st YySue, Avoop uvout 07 YySueT yyed PpeAdesqo JO O1fet OUT, 


ge e8-6IL (4) oe ey der IZl | OCF S-FF (Z “py 
oF OSE pole hy ¢-1> 
ores o1ey | : A O1yey, 
207 ¢-] Iet+ec6 (dq) L OT ‘1986 91FE-S \Z y 
(-opeo) (-opeo) (-o[eo) (-opeo) 
9/AeW UL 9/AOW UL 
nay WINYUIULOT, ee ee wn yUStOY, 
qe[d 
ai ZOT-+ / ae wu TIT 
9/AOW rer SOOT (4 “d)oV a/AOW %2~OZOT ( “)o9 (Sap) UONeZTUOT epryaed 
p OATYRIOL JO Arepuooes 
oq 0} posse AreurLg 0} VUITYSE [BUST A, 


D YOAY-oA ‘L QI, 


96 On an Unusual Double V-Event at Sea Level 


A®-particles and the results are given in tables | and 2. The errors 
quoted are standard deviations, the asymmetry occurs because the 
momentum of the primary passes through a maximum as ¢, is varied. 

It can be seen that the evidence of table 1 is consistent with particle a 
being a 0°- or a A°-particle. 

If particle 6 were a 6°-particle, the ionizations of its secondaries would 
be near to J,,,,, in disagreement with the observed ionizations. Further- 
more the secondary 6, would be visible above the plate. It is concluded 
that particle b is a A°-particle, with which view all the data are consistent. 

The angle between the planes of the V°-events is 43-8". 

As mentioned earlier the two V-particles are emitted from the inter- 
action in almost opposite directions, the lack of collinearity being 
2-54+2:°5°. 

The interpretation of the photograph is not entirely clear. ‘The 
collinearity of the V° trajectories strongly suggests the two-body decay 
of a particle at rest. However, the momentum of the centre of mass in 
the 6°+A° scheme would be 843-60 Mev/c, and in the A°+A° scheme 
826+130 mMev/c, both values obviously being inconsistent with zero. 

It is possible that the photograph shows a reaction of the type 

7-+P > P+/°, 
If the momentum of the 0°-particle exceeds 1 Gev/e at production (the 
estimated momentum of particle a is ~ 1020 Mev/c if a 0°-particle) then 
the greatest possible angle between the V-particles is 80° (James and 
Salmeron 1955). This would indicate scattering of one or both of the 
V-particles ; the low momentum, 177*{} Mev/c, of particle b is typical of 
A°-particles and may result from a loss of energy by scattering. James 
and Salmeron estimate that only 2°% of pairs of V°-particles emerge from 
a heavy nucleus without either suffering scattering. Thus double pro- 
duction with subsequent scattering seems probable, but the collinearity 
must then be fortuitous. 
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IX. CORRESPONDENCE 


_ Proportional Counter Measurements of 7-Mesonic X -Rays from Beryllium 


‘By D. West and E. F. Brapiry 
Atomic Energy Research Establishment, Harwell 


| Received October 14, 1955} 


Up to now z-mesonic x-rays have been studied exclusively with scintillation 
counters (Camac e¢ al. 1952, 1955, Stearns et al. 1954). An improvement 
in energy resolution, amounting to nearly a factor of six, can be attained 
by using proportional counters instead. An aluminium counter (diameter 
5-3 cm, wall thickness 0-030 in.) was fitted with field tubes providing a 
well-defined region, 10 cm long, over which the multiplication factor was 
constant. It was filled with xenon at 2 atm pressure, plus methane at a 
pressure of 3cm of mercury. The negative meson beam from the 
Liverpool synchrocyclotron entered the apparatus shown in fig. 1. The 
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Experimental Arrangement. A, Slowing material 1Zin. Cu; X, .. Mesonic x-ray 
source element; 1, 2, Plastic scintillators 2-3 in. x 2:5 in. x 0-25 in. ; 
3 Plastic scintillator (anticoincidence) 4 in. x 3-5in. x 0-25 in. 


thickness of the absorber A is adjusted to give a maximum number of 
7 -mesons stopped in the target material X. Scintillation counters 1, 2 
and 3 were used to select mesons stopping in X by means of a 1, 2-coincidence 
in anticoincidence with 3. The 1, 2-3 pulse was put in coincidence with 
the proportional counter output (resolving time 2psec); and this 
coincidence pulse was used to trigger a gating unit which opened for 
10ysec and allowed the proportional counter pulse to pass to a multi- 
channel kicksorter. The anticoincidence scintillator 3 completely covered 
the active region of the proportional counter and was viewed by two photo- 
multipliers. ‘The beam intensity was adjusted to give a negligible casual 
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rate. The mean rate in the proportional counter was about 1600 ¢.p.m. 
(160 000 c.p.m. during the pulse) and about 400 7~-mesons/minute stopped 
in the target. The gating rate was about 15¢.p.m. The counter was 
calibrated with K x-rays from Gd (weighted average energy of LK 
peak=41-3, kev) at intervals of two hours during the runs. 
Two runs, each of ten hours’ total duration were carried out with metallic 
beryllium samples, 2-97 g/cm? in thickness (3 in. x 2{ in.) and 2-3 g/em* 
in thickness (2} in. x 2} in.) respectively. The spectrum obtained with 
the thicker sample is shown in fig. 2. Peaks due to the mesonic K,(2p+1s) 


Fig. 2 
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a mesonic K x-rays from Beryllium (2:97 gm/cm?). 


radiation with its associated ‘escape’ peaks K,—K, yenons K-Kesecs 
are clearly seen, together with smaller peaks Kg(3p—1s) and its ‘ escape ’ 
peak Kg-K, xenon. Counts in adjacent channels were added to give the 
above spectrum. In determining the energy of the K, peak, curves were 
drawn by two independent observers using plots in which adjacent 
channels were added in pairs in the two possible ways, so that four values 
of the energy were obtained. A correction due to slight non-uniformity 
of multiplication factor along the counter is necessary since the calibration 


*K, xenon (electronic) =29:66 kev, Kg xenon (electronic) =33-62 key, 
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x-rays entered the counter mid-way along it, whereas the mesonic X-ray 
source covered most of the counter. The correction was determined by 
measuring the pulse heights from a calibration source placed at a large 
number of points over the area covered by the beryllium. It amounted 
to a change of +0:3%. Possible shifts associated with the high counting 
rate during the cyclotron pulse were checked with a strong source placed 
behind a rotating shutter. The shift was less than 0-1°% at the rate used 
during the cyclotron runs. The value for the K, energy obtained with 
the thicker source was 41-69 kev and with the thinner source 41-8, kev, 
giving a mean value of 41-7, kev. The probable error is estimated at 
+0-20 kev. Previous determinations gave values of 43-++3-5 kev (Camac 
et al. 1955), 42-9+-0-5 kev and 42:00<H<43-57 kev (Stearns et al. 1955). 
The calculated value of the energy in the absence of nuclear force effects 
is 43-9 kev (Stearns et al. 1955), so that a nuclear force shift is definitely 
confirmed for beryllium. 


Fig. 3 
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a mesonic L x-rays of Fluorine. Source LiF 1-74 g/cm?. 


The intensity of the Kg peak, which may contain contributions from 
higher transitions, is approximately 0-4 times that of the Ky peak when 
the relative efficiency of the counter for the two radiations is taken into 
account. 

The widths of the K, and its escape peak were determined by fitting 
- gaussians to the spectra obtained by plotting the single channel counts. 
Values of the observed standard deviation o,,,=1-:9 and 2-2 kev were 
obtained for the K, peak in the two runs and o,,,=1-3 and 1-1 kev for the 
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escape peak. The instrumental widths determined from measurements 
with a small source placed mid-way along the counter are 0;,,=1-3 kev 
and o;,,,=0-83 kev for these two energies. However a more direct 
measurement of the instrumental widths using the whole sensitive region 
of the counter is required. For this purpose the 7-mesonic L rays of 
fluorine were studied. Energy shifts and line broadening due to nuclear 
interactions are expected to be small in this case. A run of six hours’ 
duration was carried out with a source of lithium fluoride 1-74 g/cm? 
in thickness. The spectrum shown in fig. 3 was obtained. It clearly 
shows peaks due to the L,(32) transition and its associated escape peaks 
(La-Ky xenon) 2nd (Ly-Kg xenon). The energy was determined as 41:25 kev 
to be compared with the calculated value of 41-5 kev (Stearns et al. 1955), 
neglecting nuclear shifts but including vacuum polarization and finite 
nuclear size effects. The observed widths of the L, peak and its escape 
peak are o,,,=1:7 kev and o,,,=0-8, kev respectively, and these widths 
are presumed to be purely instrumental in origin. Since the energy is 
closely similar to that of the beryllium K radiation the fluorine spectrum 
provides good values of the relevant instrumental widths. The K, peak 
of beryllium is therefore naturally broadened by an amount 


P= VJ (Cope — Sing?) = V (2°05? — 1-72) = 1-3 kev 
whereas its escape peak is naturally broadened by an amount 
P=(1-22—0-8,2)=0-7, kev. 


Any natural broadening of the K, line would of course show up equally 
in the main and escape peaks. 

A theoretical expression for the natural broadening due to capture from 
the 1s state has been given by Brueckner (1955). In the case of beryllium 
a value of [’=0-9 kev is predicted. The indications of a broadening in the 
present measurements are consistent with this value. 
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Polygonization in Copper 


By A. Franks and D. McLran 


Communication from the National Physical Laboratory 
[Received September 29, 1955] 


RECENT experiments in this Laboratory have shown that copper tends 
to polygonize during creep. This is evident from back-reflection x-ray 
diffraction patterns taken with a microbeam camera (beam diameter 
approximately 60, specimen to film distance 14 mm). 

Figures 1, 2 and 3 (figs. 1-4, Plate 2) show the progressive development 
of polygonization in a copper specimen. The ill-defined spots in fig. 3 
indicate that the process of polygonization was not completed. Sub- 
crystal formation could not proceed beyond the stage shown in the figure 
because further strain of the specimen resulted in fracture owing to the 
brittleness of copper at high temperatures. 

Polygonization is associated with the climb of dislocations and there- 
fore the rate of polygonization depends on the number of jogs and vacancies 
(Mott 1951). The rates of polygonization in aluminium and copper were 
compared by taking microbeam photographs of a number of specimens 
which had been subjected to various rates of strain and extensions. The 
experiments were carried out at temperatures proportional to the ratio 
of the absolute melting points of aluminium and copper respectively, in 
a rough attempt to make the thermal mobilities in both metals equal. 

Figures 3 and 4 illustrate the differing rates of polygonization. The 
aluminium specimen (fig. 4) was stretched faster than the copper specimen 
(fig. 3) to approximately the same extension, but produces sharper spots. 
Though the experiments were not carried out at identical rates of strain 
and extension, the general conclusion after examining all the photographs 
is that aluminium polygonizes faster than copper. 

The difference in the rates of polygonization is however not nearly as 
great as the difference in the number of jogs present in thermal equilibrium, 
as deduced from the jog energies calculated by Seeger and Schoeck 
(1955). Using their values of about 4 ev and 0-4 ev for the energies of 
jogs in edge dislocations in copper and aluminium respectively, the ratio 
of the number of jogs present in thermal equilibrium at the temperatures 
used in the experiments (7',,= 400°C and 7’',,;= 200°C) is 

exp [(4 ev/K7'¢,,) — (0-4 ev/K7"'4;)|= 10”. 
If the number of jogs in thermal equilibrium determines the rate of poly- 
gonization, then aluminium should polygonize to the same extent as the 
copper specimen in fig. 3 in 10-1 sec. This is obviously contrary to 
experience (Wilms and Wood 1949, Greenough and Smith 1950). 

Tt can be concluded from a general examination of the present results 
that aluminium polygonizes about 100 times faster than copper. It is 
clearly unlikely that the vacancy concentrations differ sufficiently to 
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account for much of this discrepancy. The phenomenon can however be 
explained by making the assumption, which is consistent with current 
ideas, that the number of jogs produced mechanically during the deforma- 
tion of the copper specimen is of the order of 1/100 of the total number 
present in the aluminium specimen. 
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Gamma Rays from Inelastic Neutron Scattering in Sodium and 
Lodine 


By EvizaBpetu A. WOLF 
Clarendon Laboratory, Oxford 


[Received October 3, 1955] 


THE gamma ray spectrum from inelastic scattering of 2-5 Mev neutrons 
in sodium metal has been measured with a sodium iodide crystal spectro- 
meter, by the method described in a previous publication (Eliot, Hicks, 
Beghian and Halban 1954). 

The scatterer was prepared by melting small blocks of sodium under 
dry liquid paraffin, into an aluminium can 5 cm in diameter, 3 cm deep, 
and 0-1 mm in wall thickness. Any oxide that had formed rose to the 
surface of the liquid metal and could be easily removed. The container 
was then closed with an aluminium lid and sealed. To compensate for 
scattering in the aluminium an empty,can, similar to the one containing 
sodium, replaced the scatterer during background measurements. 

Pulse height distributions were measured with and without the sodium — 
scatterer, and normalized to the same number of neutrons entering the 
detector. The lower curve in fig. 1 shows the spectrum without scatterer, 
due to primary neutrons only. The upper curve shows the spectrum with 
scatterer, in which gamma rays from inelastic scattering in sodium are 
superimposed on the background due to neutrons. The spectrum due to 


sodium gamma rays only, with the background subtracted, is shown in 
fig. 2 
g. 2. 
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With fast neutrons the most likely processes to occur in a Nal crystal 
are inelastic scattering in the sodium and iodine. Gamma rays from 
scattering in sodium will be enhanced by introduction of the scatterer, 
while those from iodine will not. It can be seen in fig. 1 that several peaks 
visible in the crystal background appear stronger in the spectrum with 
the scatterer, and from fig. 2 it is evident that these are due to gamma 
rays from sodium of energy 


0-45-+0-01, 1:69-+0-03, and 2-:2+0-1 Mev. 


The remaining peaks in the background are not enhanced by introduction 
of the sodium scatterer, and are therefore most probably due to scattering 
in the iodine. They correspond to photopeaks from gamma rays of energy 
~0-2,~0-41, 0:64-40-05, and 1:01-0:05 Mev 

in good agreement with previous measurements of inelastic scattering in 
iodine (Kiehn and Goodman 1954 a, Scherrer ef al. 1954, and Lind and 
van Loef 1955). 

Previous measurements on sodium by various methods have indicated 
levels below 2-5 Mev of energy 


0-439, 2:07, and 2:37 Mev (Endt and Kluyver 1954). 


The strong 0-45 Mev gamma ray found by us is therefore almost certainly 
due to decay of the first excited state. The 1:69 Mev gamma ray may be 
due to direct decay of a new level of that energy to the ground state, or 
else to decay of a 2-14-+-0-04 mev level by cascade through the first excited 
state. The ambiguity could be most easily resolved by measuring the 
threshold energy of the 1:69 Mev gamma ray with a neutron source of 
variable energy (Kiehn and Goodman 1954). The second alternative is 
supported by the presence of the 2-2 Mev peak, which would correspond to 
an ‘ addition peak’ due to simultaneous detection of the two cascading 
gamma rays. Its intensity is somewhat greater than would be expected 
for such a peak (Kahr and Lyon 1953), but the difference is within experi- 
mental error. 
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The Interactions of Cosmic Ray Alpha-Particles 


By C. J. WADDINGTON 
H. H. Wills Physical Laboratory, University of Bristol 


[Received October 1, 1955] 


Tus letter describes the results of a study of the nuclear interactions 
produced by «-particles of the primary cosmic radiation. These particles 
were recorded in various stacks of stripped and glass-backed nuclear 
emulsions exposed at balloon altitudes. Most of the disintegrations 
were observed in a stack of stripped emulsions exposed in such conditions 
that few of the incident «-particles had energies less than 0-5 nev/nucleon : 
the median energy of all the «-particles considered was approximately 
1-3 Bev per nucleon. 

The nuclear disintegrations were detected by first making a systematic 
search for the tracks of fast «-particles (Waddington 1954), and then 
following them through the emulsions until they either left the stack, 
or interacted witha nucleus. Highty-nine ‘stars’ were observed in a total 
effective path length of 1825cm; the corresponding value of the inter- 
action mean free path of an «-particle in emulsion is 20-5+2-2 em. 
A separation of the tracks into groups corresponding to «-particles of 
different energies, failed to reveal any significant variation over the energy 
range considered. Further, the mean value is in good agreement with 
that found by Quareni and Zorn (1955) for the mean free path of 350 Mev 
a-particles, i.e. 20:02+1:67 cm. In only two examples was a high energy 
a-particle observed to emerge from one of the disintegrations, a result 
indicating that the mean free path for absorption is only a little greater 
than that for interaction ; if the «-particle interacts with a nucleus it 
rarely emerges. 

The interaction cross section corresponding to the observed value of 
the mean free path, as observed in the present experiment, is equal to 
approximately half the value calculated from the equation derived from 
geometrical considerations :— 


et el), 


where R, and R, are the radii of the target and incident nuclei 
respectively, their magnitudes being given, for a nucleus of atomic 
number A, by the relation R=r,A¥%, with rp=1-45x10-% em. If 19 is 
taken as 1-20x10-!%cem, the disparity between the observed and 
calculated values is still considerable: o,,=0-70c,,. If, alternatively, 
the cross section is calculated assuming R,=0, it becomes greater than 
the experimental value, even when 7, is taken to be 1-45x 10-** cm, 
Since the experimental cross sections found for energetic nucleons, and 
for the heavy nuclei of the cosmic radiation, are always less than those 
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derived geometrically, the assumption that R, may be neglected does 
not seem justified. ina 
Sixty-eight of the stars observed in the present investigation, five 
of which were tentatively identified as due to collisions with hydrogen, 
have been analysed, in order to compare them with those observed by 
Ceccarelli ct al. (1955). These authors detected the disintegrations 
due to «-particles by a direct search. It is well known that this method 
commonly suffers from the fact that there may be a considerable difference 
in the efficiency with which observers detect disintegrations of different 
types. In particular, there is a tendency to overlook small stars, a 
defect from which the present method is free. A comparison between 
the results from the two methods allows the effect of any scanning 
bias to be estimated. Figure 1 shows the distribution in the values 
of N, for the stars observed in the two experiments. It has been 
convenient to plot an integral distribution—the number of ‘stars’ 
with NV, greater than any given value—and to normalize the two values 
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The integral NV, distribution for the stars found in the present investigation 
is shown by the continuous line; and that found by Ceccarelli e¢ al. 
by the broken line, the two distributions being normalized at N,—=7. 


corresponding to N,—7. The results suggest that there was a consi- 
derable bias against detecting small stars in the earlier work. Whereas 
the early data suggested that only about 2°, of the ‘stars’ had N aw‘, 
the present results indicate that the true value is about 50% (see also 
Engler et al. 1954). 

A similar indication of the effects due to a failure to observe small stars, 
in methods employing direct scanning, is provided by the results shown in 
fig. 2, which gives the distribution in n, of the stars found in this present 
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investigation. The results of Ceccarelli et al. indicate that there is 
a peak in such a distribution at n,=1, whereas in the present work 
it is actually observed at n,=2. This discrepancy can be readily 
accounted for if the majority of the small stars with N,<4—believed, 
on the basis of the evidence presented in fig. 1, to have been overlooked 
in the early work, have two associated shower particles. This interpreta- 
tion is supported by the distribution in the values of n, of those stars 
observed in the present experiments which show a total liberation of 
energy, among the evaporation and ‘ knock-on’ fragments, of less than 
100 Mev, and which have not been tentatively classified as collisions 
with hydrogen. The values of n, for such stars, all of which have 
N,<4, are shown shaded in fig. 2, and display a very pronounced peak 
at n,=2 which accounts for most of the peak observed in the total 
distribution. 


Fig. 2 
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The distribution in ”, of the stars found in the present investigation. Those 
shown ‘ blocked-in’ correspond to ‘stars’ which, while showing less 
than 100 Mev visible energy in the form of black or grey ‘ prongs’, 
have not been classified as due to collisions with hydrogen nuclei. 


The observed features—both of the mean free path and of the 
distribution in N, of the disintegrations—may be explained by dividing 
the interactions into two separate classes. For simplicity, this may 
be done on the basis of a purely geometrical model, in which the significant 
parameters are taken to be the overlapping target areas of nuclear 
matter in the collisions rather than the interpenetrating volumes. The 
two types of collisions are then : 

(i) Central collisions in which all the nucleons of the incident «-particle 
pass through the nuclear matter of the target nucleus. The interaction 
cross section for such collisions will be given by 

g,=7(h,—f,)*. 
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(ii) Peripheral collisions in which not all the nucleons of the incident 
nucleus pass through the nuclear matter of the target nucleus. 

The relative cross sections for collisions of these two types is independent 
of the assumed value of 7, and it is found that for emulsion nuclei only 
~14°%, of the total cross section will be available to central collisions. 
The results of Ceccarelli et al. indicate that for these central collisions 
the nucleus is opaque, even when allowance is made for the transparency 
of the nuclei to individual nucleons. In order, therefore, to account 
for the small value of the experimental cross section, it is necessary to 
assume that the effectiveness of peripheral collisions is small. For 
such collisions, the nucleus is partially transparent. The degree of 
transparency deduced from the results depends on the assumed value of 
79, and represents an average over the various possible impact parameters. 
Using the experimentally determined cross section, and assuming that 
p= 1-45 x 10-8 em, the transparency f, for peripheral collisions is found 
to be 0-43 ; the corresponding value for the fraction of central collisions 
among those observed is 27%. If, however, 7>=1-20x 10~¥ cm, then 
the transparency is f=0-71, and the percentage of central collisions is 
18°, of those observed. 

In the present work, the number of stars with N,,>9 was 27% of the 
total. Since the cross section for central collisions is almost entirely 
due to the heavy nuclei of the emulsion, it appears reasonable to assume 
that the number of such stars represents an upper limit to the number 
of central collisions. This assumption is supported by the observed 
peak at n,=2 in the distribution shown in fig. 2, and at n,=1 in the 
earlier work in which, it is suggested, there was a much larger proportion 
of central collisions. 

It seems therefore that the simple model of a ‘ black ’ nucleus proposed 
by Ceccarelli e¢ al. should be replaced by one in which the core only of the 
nucleus is opaque, while the edges have a mean transparency of which 
the estimated value depends on the value taken for 7). The use of such 
a model, in calculating the mean free paths in media other than emulsion, 
gives results which are in good agreement with the experimental values 
in glass, found by Bradt and Peters (1950); and with the values based 
on their model, which assumes that a certain overlap of the nuclear radii 
must oceur before an interaction can take place. 


REFERENCES 
Brabr and Prrurs, 1950, Phys. Rev., 77, 54. 
CECCARELLI, QUARENI and Zorn, 1955, JI Nuovo ¢ “m., 1, 669. 
Enouer, HArBEeR-ScHaim and WINKLER, 1954, J1 Nuovo Cim., 12, 930. 
QUARENI and Zorn, 1955, Il Nuovo Cim.. 1, 1282. 
WappvinaTon, 1954, Phil. Magq., 45, 1312. 


Correspondence 109 
The Paramagnetic Resonance Spectrum of Chlorine Dioxide 


By J. KE. Bennerr and D. J. E. Ineram 
University of Southampton 


[Received October 26, 1955] 


CHLORINE DIOXIDE is one of the few stable gases which possess unpaired 
electrons in their normal chemical state and are therefore paramagnetic. 
Detailed resonance measurements have been made on the other common 
gases, such as nitric oxide, nitrogen peroxide and oxygen (Beringer and 
Castle 1950, 1951), but none have so far been reported on chlorine dioxide. 
Preliminary experiments (Bennett e¢ al. 1955) have indicated, however, 
that an intense absorption could be obtained from aqueous solutions 
of chlorine dioxide with a marked hyperfine structure from the nuclear 
moment of the chlorine atoms. 

In this letter measurements on different solutions of chlorine dioxide 
are reported and the variation of the hyperfine splitting with the physical 
state of the solvent shows that a strong ‘averaging’ or ‘ tumbling ’ 
effect is produced by the mobile solvent molecules. The simplest type 
of spectra is that obtained from a dilute solution of chlorine dioxide in 
alcohol at 1-25cm wavelength and 200°K, as shown in fig. 1. Four 
narrow and well resolved lines are observed, the width of each line being 
8 gauss with an overall separation of 50 gauss. The g value of 
2-010-+-0-001, corresponding to the centre of the pattern, can be measured 
by comparison with the hydrazyl absorption line, which is exactly 

‘coincident with the highest field hyperfine component at 1:25 ecm 
wavelength. When this same solution is cooled to liquid oxygen tem- 
peratures a marked change in the spectrum is observed; as shown in 
fig. 2. It would appear that the hyperfine splitting has considerably 
increased and also become markedly anisotropic. As a result the 
individual lines are noticeably broadened and the pattern corresponds 
to an averaging over all possible orientations, with a maximum splitting 
of 170 gauss, which is more than three times that observed in the liquid 
solvent. 

These results may be most readily interpreted by assuming that there 
is a strong crystalline field axis acting in each molecule when in the solid 
solution, which produces the anisotropy in the hyperfine splitting. 
The very narrow and well resolved lines obtained from the liquid solution 
show that the motion of the solvent molecules is extremely effective in 
averaging out the anisotropy in the splitting. This radical change in 
the spectrum on freezing the solvent is quite general as the two different 
patterns are also obtained from solutions in acetone, benzene and. water. 
The temperature at which the pattern changes always corresponds to 
the melting point of the particular solvent. 

The lines in the liquid solvent are broadened by both * spin-spin ” 
and ‘ spin-lattice ’ interaction, Thus, if the concentration of the dissolved 
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chlorine dioxide is too high, the individual hyperfine components are 
not resolved at any temperature and a single broad absorption of 150 gauss 
total half-width is produced. This same line is also obtained from the 
pure chlorine dioxide in either the liquid or solid phase. A broadening 


Fig. 1 


ee 


Spectrum observed from alcoholic solution of chlorine dioxide at 200°K and 
1-25 cm wavelength. 


Fig. 2 


——— 


Spectrum observed from same solution as in fig. 1, but at 90°K. 


The two spectra are on the same scale and the line at the bottom is equivalent 
to a 100 gauss splitting. 


of the individual lines is also observed as the temperature is raised, 
and in the case of alcohol only a single unresolved line is obtained at 
room temperature. This, together with the departure of the g value 
from that of a free-spin indicates that there is considerable orbital 
coupling to the unpaired electron, 
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Another feature of interest is the difference between the hyperfine 
patterns at 38cm and 1-:25cm wavelengths. The solution which gives 
the four lines of fig. 1 at 1-25 em wavelength gives an identical pattern 
at 8mm, but only a single line of 26 gauss width is obtained at 3 cm 
wavelength. This is again true for all solvents, and it would appear 
that a form of exchange narrowing is present so that frequencies higher 
than 10000 Mc/s are necessary to resolve the individual hyperfine 
components. In. contrast to this the same spectra is observed at all 
wavelengths at 90°K showing that the ‘tumbling’ effect of the mobile 
solvent molecules, and any exchange interaction between the chlorine 
dioxide molecules, have both been quenched at this temperature. 
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X. REVIEWS OF BOOKS 


Radio Astronomy. By J. L. Pawsry and R. N. BRacEWELL. (Oxford : 
University Press, 1955.) [Pp. x+357+23 plates.] Price 55s. 


To present a well-balanced review of a subject which is expanding as rapidly 
as radio astronomy is no easy task and the authors of this volume, another 
in the series of International Monographs on Radio, have achieved a commend- 
able degree of success. In this book, which is certainly the most comprehensive 
account of the subject which has yet appeared, the growth of the various 
branches is taken to approximately mid-1952. It was during this period that 
most of the fundamental advances were made; these include the accurate 
location of the most intense radio stars and their subsequent identification 
with visual objects both inside and outside the galaxy, the discovery of the 
21 cm hydrogen line, the detailed studies of thermal radiation from the solar 
corona and the spectrum analysis of the enhanced radiation associated with 
solar activity. 

The book appears to be written primarily for the physicist with no previous 
knowledge of the subject, and the first three chapters are devoted to an account 
of the observational techniques and to the theory of radio waves in ionized 
gases. Two additional chapters deal with some of the basic features of general 
astronomy and astrophysics which are relevant to radio astronomy. Further 
chapters give a detailed account of solar, lunar and cosmic radio emission, the 
investigation of meteors and the moon by the radar echo technique, and effects 
introduced by the earth’s atmosphere. Throughout the book the reader is 
constantly referred to original papers and detailed references are given at the 
end of each chapter. The text is well illustrated by line diagrams and some 
two dozen plates are included. The latter show different types of aerial system 
and also some photographs taken at the Mt. Wilson and Palomar observatories. 
Unfortunately we are not shown the colliding galaxies in Cygnus. 

Impeccably printed by the Oxford University Press this book can be warmly 
recommended to all who desire a general survey of the principal branches of 
radio astronomy, and it should be of particular value to research workers who 
are newcomers to the subject. A. H. 


Applied Group-Theoretic and Matrix Methods. By B. Hiteman. (Oxford: 
Clarendon Press.) [Pp. xii+454.] Price 60s. 


TxIs book gives an account of the applications of group theory to a wide 
variety of subjects in physics and chemistry including some subjects less 
frequently treated such as Eddington’s quantum relativity. It is useful as 
giving an indication of the fields in which the theory can be applied, but in 
covering so wide a field it seems to the reviewer to fall between two stools : 
the treatment is not rigorous and not complete enough to serve as an authorita- 
tive work on any one branch of the subject ; it is on the other hand too advanced 
and too long to serve as an introduction to the subject for a beginner who wants 
to learn it with a view to its applications in natural science. Also the book 
makes use of a notation for vectors which is different from that most commonly 
in use, and the explanation of the crystalline notation is spread over a very 
considerable part of the book. The presence of many small defects such as 
these would make it difficult to use as a reference book, Veo 
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Diffraction photograph of annealed copper. 

Diffraction photograph of copper stretched 6—7°% in 57 hr at 400°c. 
Diffraction photograph of copper stretched 6°, in 308 hr at 400°C. 
Diffraction photograph of aluminium stretched 8-4°% in 5 hr at 200°c, 


